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study of Radiative Decays of Psi(2S) Mesons 

by Kaiyan Gao 
Under the supervision of Professor Ronald Poling 

ABSTRACT 

We studied the decay ^(2-5) 7 r]c{2S) with 25.9 x 10^ ^(2-5) events collected 
with the CLEO-c detector. No il^{2S) — > 7 ric{2S) decays were observed in any of 
the eleven exclusive r]c{2S) decay modes studied. The product branching fraction 
upper limits were determined for all modes. The 90% confidence level upper limit of 
B{i;{2S) 7 Vc{2S)) < 7.4 X 10"^ was obtained. 
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CHAPTER 1 



Introduction 

This dissertation is devoted to the study of the radiative decay of charmonium me- 
son ip{2S) created in e~^e~ annihilations at the center of mass energy of 3.7 GeV. 
Specifically, we have studied the radiative decay ip{2S) '-yrjci^S) using CLEO-c 
detector by reconstructing the final states of r]c{2S) decays exclusively. To provide 
the background for the research I start with the fundamental knowledge of particle 
physics. 

Since the start of human history, the composition of the world has always been 
a subject of intense interest. From thousands of years ago, when the name atom 
was given, to recent decades, experimental observations and theoretical studies have 
allowed us to understand the fundamental nature of materials more and more deeply. 
The latest and the most successful model of the fundamental particles that compose 
the universe is the Standard Model, which began with the introduction of quarks by 
Murray Gell-Mann and George Zweig in 1964. The observation of the top quark 
completed the quark sector of the Standard Model [H [2] and provided strong support 
for the model. It is not a complete description, however, as demonstrated by recent 
observations of new physics beyond the Standard Model, which may greatly expand 
our knowledge about the universe. 
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1.1 The Standard Model 



Nowadays it is common knowledge that all materials are composed of atoms and 
molecules. Atoms are formed from elementary particles: electrons, protons and 
neutrons. Since the discovery of the first elementary particle, the electron, in 1897 
by J. J. Thomson at the Cavendish Laboratory of Cambridge University, hundreds of 
elementary particles have been found and various theoretical models have been built 
to explain the properties and interactions of these particles. For over forty years, the 
Standard Model has been developed as a framework for understanding and studying 
these elementary particles. Experimental measurements have demonstrated that 
the Standard Model provides an excellent description and have provided increasingly 
precise determinations of many unpredicted parameters. 

1.1.1 Quarks and Leptons 

Within the Standard Model, the particles with no known substructure are categorized 
in three types: quarks, leptons and gauge bosons. The quarks and leptons are the 
building blocks of composite particles, and are fermions with spin s — ^. The gauge 
bosons are the mediators of the fundamental interactions and have whole- integer spin. 

Hadrons are composite particles made of quarks. There are six kinds of quarks 
and they are characterized according to their "flavor": u, d, c, s, t, and b (up, down, 
charm, strange, top, and bottom quarks) with charges (in units of the electron charge) 
of |, — |, |, — |, |, and — |, respectively. The flavor quantum numbers are the third 
component of isospin of u and d quarks, the charm of c quarks, the strangeness of 
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s quarks, the topness of t quarks, and the bottomness of b quarks. 

The properties of the quarks are summarized in Table ll.li Quarks have not 
been found to exist independently and are only found as constituents of composite 
particles. The explanation of this is that quarks have another quantum number, 
color, so named because of similarity with color theory of visible light. Because net 
color is not observed in nature, the quarks are "confined" into the composite particles 
called hadrons, which are colorless. 

Table 1.1: Properties and additive quantum numbers of the quarks: mass (M), 
electric charge (Q), isospin 3-component {I3), Charm (C), strangeness (5), topness 
(r) and bottomness(i?). Three generations of quarks are separated by thick lines. 



Flavor 


M(GeV/c2) 


Q 


h 


c 


s 


T 


B 


u - up 


3 MeV/c^ 


^3 


+\ 














d - down 


6 MeV/c2 


1 
3 


1 
2 














c - charm 


1.24 GeV/c2 


^3 





+1 











s - strange 


95 MeV/c2 


1 
3 








-1 








t - top 


172 GeV/c^ 


^3 











+1 





b - bottom 


4.2 GeV/c^ 


1 
3 














-1 



There are two kinds of hadrons. Baryons consist of three quarks and are therefore 
fermions, e.g. a proton p ~ uud and a neutron p ~ ddu. Mesons consist of one quark 
and one antiquark and are therefore bosons, e.g. 7r+ ~ ud and ~ us. Hadrons 
are particles that participate in strong interactions directly. 
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The Standard Model leptons are the electron (e~), muon and tau (r^), each 
with charge of —1, and their corresponding chargeless neutrinos [ug, and v^-). Some 
of the properties of the leptons are given in Table fL2] Each lepton type (e~, yU~, r^) 
is associated with a conserved lepton number. Recent evidence of oscillations among 
neutrino flavors and nonzero neutrino masses provides the first evidence for physics 
beyond the Standard Model P H El E] . 

Table 1.2: Properties of the leptons. Three generations of leptons are separated by 
thick lines. 



Flavor 


Mass 


Electric Charge 


z/g - electron neutrino 


< 2 eV/c^ 





e - electron 


0.511 MeV/c^ 


-1 


- muon neutrino 


< 0.19 MeV/c2 





^ - muon 


106 MeV/c^ 


-1 


Vr - tau neutrino 


< 18.2 MeV/c^ 





r - tau 


1.78 GeV/c2 


-1 



Tables 11.11 and 11.21 show the natural grouping of the quarks and leptons in three 
generations. Why there are three generations is one of the deepest mysteries of 
particle physics. 

For each quark and lepton there is a corresponding antiquark and antilepton, and 
as described above the quarks come in three colors. The physical properties of the 
quarks and leptons from different generations are identical except for the masses of 
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the fermions. 

1.1.2 Intermediate Particles and Fundamental Interactions 

The Standard Model describes the electromagnetic, weak and strong interactions. It 
does not incorporate the fourth of the fundamental interactions, gravitation. Each in- 
teraction has its corresponding "force carrier" , particles that are called gauge bosons. 
These mediators are all bosons with spin s — 1. Massless photons (7) are chargeless 
and associated with the electromagnetic interaction, which has infinite interaction 
range. Massive PF^'s and Z°'s are associated with the weak interaction, which has a 
very short range. They are self- interacting and the W^''s have electric charge of ±1, 
while Z°'s are neutral. Eight massless and chargeless gluons (g) are associated with 
the strong interaction, which is also short-ranged because the gluons carry color. As 
gluons have eight different colors, they not only interact with quarks but also interact 
with themselves. 

1.1.3 Parameters of the Standard Model 

A fundamental feature of the weak interaction in the Standard Model is quark mixing. 
The weak eigenstates of quarks are not the same as their mass eigenstates. The 
mixing of flavors and generations induced by weak interactions is parameterized by 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, V, a 3x3 unitary matrix, operating 
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on quark mass eigenstates, as shown in Equation ll.li 





d' 


\ 




Vud 


Vus 


Vub 


\ 




d 


\ 






s' 






Vcd 




Vch 






s 




(1.1 


v 


h' 


I 


\ 


Vtd 


Vts 




1 


\ 


b 


I 





The 2008 Nobel Prize in physics was awarded to M. Kobayashi and T. Maskawa for 
this formulation, "the discovery of the origin of the broken symmetry which predicts 
the existence of at least three families of quarks in nature". Without flavor mixing, 
quarks could only decay within their generation. We would expect the heavier 
member to decay 100% of the time to its lighter partner (c ^ s, t ^ 6) and the 
lighter one to be stable. Since the CKM matrix is unitary, the nine elements are 
reducible to four independent parameters. In the standard parameterization of the 
matrix, three angles (6'i2, 6'23, Oi^) and one phase (^13) represent the four independent 
parameters. 

In experimental high energy physics, in addition to the properties of particles such 
as masses, charges and spins, the fundamental parameters of the CKM matrix must 
also be measured. These parameters are reflected in the decay widths (lifetimes), 
branching fractions, and other detailed properties of the particles and their decays. 
Redundant measurements of the CKM parameters are important and necessary as 
they are a powerful test on the viability of the Standard Model. 

Based on various measurements, the current ranges for the magnitude of the CKM 
matrix elements have been compiled by the Particle Data Group [7] and are given in 
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Equation II. 2t 



0.9741 



0.9756 



0.219 



0.226 



0.0025 



0.0048 



\ 



V 



0.219 



0.226 0.9732 



0.9748 



0.038 



0.044 



(1.2) 




0.014 0.037 ^ 0.044 0.9990 ^ 0.9993 ) 



The fact that the couphng between different generations of the quarks is small can 
be seen from the small values of the non-diagonal elements. 

1.1.4 Symmetries and Conservation Laws 

If a system or phenomenon remains unchanged under a transformation of one or more 
physics variables, then the system or the phenomenon is regarded to have symmetry 
with respect to this transformation. According to Noether's Theorem [8], any sym- 
metry of a physical system under a transformation that does not explicitly depend on 
time has a corresponding conservation law. The interactions and decays of particles 
are governed by conservation laws. 

One type of conservation law is called an exact conservation law. The most 
commonly used conservation laws, conservation of energy, conservation of linear mo- 
mentum, conservation of angular momentum, and conservation of electric charge are 
all exact conservation laws. Conservation of energy is associated with symmetry 
under time translation. Conservation of linear momentum is a mathematical conse- 
quence of symmetry under continuous translation in space. Conservation of angular 
momentum corresponds to the continuous rotational symmetry of space. 

Another type of conservation law is an approximate conservation law. Such a law 
holds only under some particular situations. 
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Isospin 

Isospin symmetry is an important internal symmetry that was introduced early in the 
development of nuclear and particle physics. The concept of isospin came from the 
similarity between neutrons and protons. Neutrons and protons have the same spin 
of I and nearly the same mass, but neutrons are chargeless and protons have +1 unit 
of electric charge. They behave nearly identically under the strong interactions but 
differently in weak and electromagnetic interactions. By introducing isospin space, 
the abstract internal space associated with isospin quantum number /, neutrons and 
protons can be regarded as a doublet of the quantum state with I — \- The algebra 
of isospin is the same as that of angular momentum and isospin is also an additive 
quantum number. Therefore protons and neutrons have different projections on 
the third dimension in isospin space, i.e. 1^ — |, — | for the proton and neutron, 
respectively. Similarly, pions have isospin of 1, so by following the algebra of angular 
momenta, they have three possible third components of isospin, — —1,0,1 and 
the corresponding pion states are denoted as tt", 7r°, tt"^. The third component of 
isospin is chosen by convention based on the charge of the particle, e.g. the proton 
has charge +1 and is assigned the positive value of I2,. 

Isospin symmetry is an approximate symmetry describing the strong interaction. 
The conservation of isospin requires that a system remains unchanged in strong in- 
teractions, i.e. the isospin / and its third component are conserved. Similar 
to the relationship between the conservation of angular momenta and the rotational 
symmetry of space, the conservation of isospin comes from the rotational invariance 
of strong interactions in isospin space. Conservation laws of isospin does not hold in 
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weak interactions. 



Parity 

Parity (V) is a discrete (non-continuous) operation and is defined as the inversion of 
space coordinates: 

X ^ -X. (1.3) 

According to its definition, parity satisfies V'^l'i/j] — ■0- Therefore, it has two eigen- 
values, P — ±1. Under P, the signs of both position x and hnear momentum p 
are reversed, but angular momentum, defined asL = xxp = L remains unchanged. 
Thus X and p have parity of —1 and L has parity of The classical variables, the 
time, the energy, the angular momentum of a particle, and masses, charges, coupling 
constants, etc. satisfy P — +1. The variables with P — —1 include the position, the 
velocity, the linear momentum of a particle, and electric field. 

The invariance under the V transformation is associated with conservation of 
the parity quantum number. Parity is conserved in strong and electromagnetic 
interactions. 

Parity is a multiplicative quantum number, i.e. the parity of a multiparticle state is 
the product of the parities of all component particles. The wave function of a particle 
is characterized by its intrinsic parity, which is the parity quantum number P of the 
particle. Basic nonrelativistic quantum mechanics demonstrates that if the orbital 
angular momentum of a system is L, then the parity of the system is the product 
of (—1)^ and the intrinsic parity. The intrinsic parities of particles are defined in 
accordance with the conservation law, and are assigned based on experiments and by 
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convention. By convention, quarks have positive parity and antiquarks have negative 
parity. The intrinsic parity of a hadron can be determined by the strong interactions 
that produces it, or through decays not involving the weak interaction. The intrinsic 
parity of a neutral particle can be determined from experiments, e.g. P(7r°) = — 1. 
Photons and gluons have ^(7) = —1 and P{g) = —1, respectively. The intrinsic 
parity of a particle composed of a fermion and antifermion pair is —1 and the intrinsic 
parity of a particle composed of a boson and antiboson pair is +1. 

Chcirge Conjugation 

Another discrete transformation is the charge conjugation (C) transformation, which 
turns a particle into its antiparticle. The charge conjugation operator C reverses 
the electric charge and all the internal quantum numbers, including lepton number, 
baryon number, /a, strangeness, and the other flavors. The eigenvalue oiC, C — ±1, is 
called charge conjugation parity (C), similar to P parity. Only truly neutral systems 
with all quantum charges and magnetic moment zero are eigenstates of charge parity. 
Photons and gluons both have C parity —1. For a system that is composed of a 
particle and antiparticle pair, if the orbital angular momentum is L and the sum of 
the spin angular momenta of the two particles is S, the C parity is C = (—1)^+'^. C 
parity is also a multiplicative quantum number. 

Invariance under the charge conjugation transformation leads to the conserva- 
tion of charge parity and is called C-symmetry. C-symmetry holds in strong and 
electromagnetic interactions, and even gravitational interactions. However, weak 
interaction does not obey charge conjugation symmetry. 
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Other Conservation Laws 

There are other conservation laws, e.g. conservation of strangeness and baryon num- 
ber, that also govern elementary particle interactions. The strangeness S, defined 
as the number of strange antiquarks s minus the number of strange quarks s, is a 
quantum number introduced to describe decays of particles in strong, electromag- 
netic and weak interactions. It has been observed experimentally that the change 
of the strangeness quantum number satisfies AS = for strong and electromagnetic 
interactions and \AS\ = 0, 1 for weak interactions. 

Baryon number is defined as one third of the difference between the number of 
quarks q and the number of antiquarks q, where the fraction one third comes from 
the fact that a baryon consists of three quarks or antiquarks. The conservation of 
baryon number holds in all interactions of the Standard Model. Lepton number is 
similarly conserved in all interactions. 

More detailed descriptions and discussions about conservation laws can be found 
in any standard textbook, such as Ref. [9]. 

1.2 Charmonium States 

Quarkonium is a fiavorless bound state which consists of a quark and its own antiquark 
bound by strong interactions. Quarkonia include charmonia and bottomonia. A 
charmonium meson(cc meson) is the bound state of a charm quark and a charm 
antiquark, and a bottomonium meson {bb meson) is the bound state of a bottom 
quark and a bottom antiquark. No toponium exists because of the large mass of 
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top quark and antiquark. The quark would decay through electroweak interactions 
before a bound state can form. Bound states of hght quark and antiquark pairs 
are normally not called quarkonium, partly because in the experiments the states 
observed are usually the mixture of the light quark states. The large mass difference 
between the charm or bottom quark and the lighter quarks prevents the charmonium 
and bottomonium states from mixing with each other or with the lighter flavorless 
mesons. Charmonium and bottomonium are very important for studying the strong 
interaction. The two charmonia involved in this dissertation are the ip{2S) (often 
written as ip' or '?/'(3686)) and rjc{2S) (often written as rj'^. 

The charmonium resonances were first discovered in experiments in 1974 [TOl [TTl 
[12], known as the "November Revolution," when the J /if) meson was independently 
discovered by a group at the Stanford Linear Accelerator Center and by a group at 
Brookhaven National Laboratory. The charmonium resonances were quickly inter- 
preted as cc bound states and were very powerful in establishing the reality of quarks 
and the Standard Model. 

Non-perturbative Quantum chromodynamics (QCD) is used to compute the prop- 
erties of mesons. Similar to nonrelativistic models of the hydrogen atom, the motion 
of the quarks in a quarkonium state is nonrelativistic so they can be assumed to move 
in a static potential. One technique that can be effectively applied to quarkonium is 
application of quantum mechanics, with specific choices of the quark- antiquark po- 
tential. Such calculations predict the masses of quarkonia states, and can be directly 
tested with data. One of the candidate potentials is [13] 

^(r) = -- + ^, (1.4) 
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where r is the effective radius of the quarkonium state and n and a are parame- 
ters. The first term is the Coulomb-type force that dominates at the short distances 
(asymptotic freedom) . It corresponds to the potential induced by one-gluon exchange 
between the quark and its anti-quark, analogous to the Coulomb electromagnetic po- 
tential. The second term is linear in r and produces quark confinement. 

The charmonium states, as well as other elementary particles, are denoted with 
the spectroscopic notation n^^^^Lj. The principal quantum number n = 1,2,3... 
represents the ground state, and successive excited states of the particle. S represents 
the spin of the particle, which for charmonium takes values S — 0,1, since its quark 
and antiquark constituent are fermions and have spin projections of ^ or — i. L 
denotes the orbital angular momentum and the letters S, P, D ... stand for L — 0,1, 2, 
... . J is the total angular momentum with J = |L -|- S"!, |L -|- 5*1 — 1, \L — S\. 
Specific states are frequently denoted by their J^^ , where P and C are the parity and 
charge conjugation quantum numbers, respectively. For example, the two charmonia 
involved in this dissertation are r]c{2S) and ip{2S). The 25" means that for these 
states, n — 2 and L — 0. In fact, rjc{2S) denotes l-'^5'o and '^{2S) is I^Sq, so rjc{2S) 
has S = and ip{2S) has S = 1. The total angular momenta are J = for r]c{2S) 
and J = 1 for ■0(25'). According to the parity convention for quarks and antiquarks, 
both r)c{2S) and ^(25*) have parity —1. However, the C parity of rjc{2S) is -|-1 and 
that of '4>{2S) is —1 because C — (— l)^"*"*^. Thus the characterization of rjc{2S) is 
jpc ^ Q-+ ^hat of ^{2S) is J^^ = 1—. 

The charmonium state r]c{2S) and ip{2S) are the first excited states of r]c{lS) 
and J/ 1^(13) respectively. Therefore the transition between ric{2S) and ■0(25) is 
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similar to the transition between ric{lS) and J/ip{lS). The lowest energy states for 
mesons with quark- ant iquark pairs have S* = and negative parity and are called 
pseudoscalar mesons, e.g. rjci^S) and rjc{2S). For excited states of mesons the quark 
spins are parallel, so if a state has zero orbital angular momentum, then the total 
angular momentum of the particle is 1. Such states are called vector mesons, e.g. 
J/ip and iIj{2S). As the vector mesons have the same quantum numbers as photons 
(virtual photons), i.e. J^'" = 1 , they can be produced directly in e~^e~ collisions, in 
which virtual photons are formed. Other charmonium states can be observed in the 
transitions of J^'" = 1 states or in pp collisions. Figure 11.11 shows the spectrum 
of the charmonium resonances and the observed transitions among them. 



1.3 r]c{2S) Production at ip{2S) 



1.3.1 ij{2S) Decays 

The charmonium states with J^*" = 1 have four possible ways to decay, as shown 
in Figure II. 2[ In experiments, they are observed to be through leptonic decays, 
hadronic decays, and radiative decays. 



(a) Leptonic decays (Figure 1.2(a)): The quark and antiquark pair annihilate to 
create a virtual photon and then the virtual photon produces a lepton and 
antilepton pair, i.e. J/ilj,ip{2S) — > 7* 



(b) Hadronic decays by the electromagnetic interaction (Figure 1.2(b) ): The quark 



and antiquark pair annihilate to create a virtual photon and then the virtual 
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0140406-002 




Figure 1.1: The spectrum of charmonium resonances. The states are labeled with 
spectroscopic notation n'^^'^^Lj, where n is the principal quantum number, S rep- 
resents the spin of the particle and S = 0,1, L = S,P,D,... denotes the orbital 
angular momentum of L = 0,1,2, and J is the total angular momentum. In 
addition to this notation, parity (P) and charge conjugation (C) are used in the 
notation J^^ . 
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photon produces a quark and antiquark pair. The quark-antiquark pair frag- 
ments into a final state of hadrons, i.e. J /ip., ^/'(2S') ^ 7* — > hadrons. 



(c) Hadronic decays by the strong interaction (Figure 1.2(d)): The quark and 
antiquark pair annihilate to three gluons, then the three gluons fragment into 
a final state of hadrons, i.e. J /ip, ipi'^S) ^ g + g + g ~* hadrons. 



(d) Radiative decays (Figure 1.2(c)): The quark and antiquark pair annihilate to 
a photon and two virtual gluons, then the two gluons produce the final state of 
hadrons, i.e. J/ip, ipi^S) ^7 + 5' + 5'^7 + hadrons. 

The strong interaction dominates over the electromagnetic interaction, in spite of 
the suppression (OZI) of the three-gluon process. It is impossible for a charmonium 
state with J^*" = 1 to decay to a single gluon as the particle is colorless but a gluon 
carries color charge. Decay through two gluons is forbidden by C conservation. For 
example, 97.85% of 'ip{2S) mesons decay to hadrons through three gluons and only 
1.73% of ip{2S) mesons decays to hadrons through virtual photons. 



1.3.2 Radiative Decay of t/j{2S) jr]c{2S) 

The meson ip{2S) has a mass of (3684 ± 0.034) MeV and can decay radiatively to 
another cc meson, such as ipi^S) -yrjci^lS) and ip{2S) lXc2- The radiative decay 
ip{2S) — > 7?7c(2S'), that is the focus of this dissertation, has not been observed. 

In our investigations, ip{2S) mesons are produced by the annihilation of electron- 
positron pairs at a center-of-mass energy of 3686 MeV. An unknown portion of the 
^^{28) mesons would decay into a transition photon and a pseudoscalar cc meson 
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(a) V(25) ^ 7 ^ l+l- (b) V(2-S) ^ 7 ^ hadrans 

7 




Hadrons 



C C 



(c) ip{2S) g + g + g ^ hadrons (d) ip{2S) -^j + g + g^^ + hadrons 

Figure 1.2: cc annihilation and decay mechanisms for tp{2S) {ip'): (a) leptonic decay 
through virtual photon; (2) hadronic decay through virtual photon; (c) hadronic 
decay through three gluons; (d) radiative decay through two gluons and a photon. 
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'rjc{2S) through electromagnetic interaction. Some of the generated 'rjc{2S) mesons 
would decay into hadrons, ultimately into detectable particles, like pions, kaons, and 
photons. The first step of the radiative decay is a two body decay, which means the 
particle rjc{2S) could be observed and measured through the energy distribution of 
the transition photon. 

The experimental challenge of the measurement for this decay channel is to detect 
the ~ 48 MeV radiative photons in an experimental environment with considerable 
background. Inclusive study is impractical for these low photon energies. Our 
alternative approach is to study the decay exclusively by reconstructing rjci^S) from 
its hadronic decay products and looking at the r)c{2S) candidate mass for evidence of 
iP{2S)^^ri,(2S). 

1.3.3 Theoretical Background of ^^(25) 7 r]ci2S) 

The purpose of this analysis is to study the direct magnetic dipole (Ml) radiative 
transition '4>{2S) — > 7 r]c{2S) by reconstructing exclusive decays of the r]c{2S). The 
partial width for a direct Ml radiative transition between S-wave charmonium states 
is given by 

res, 7) = 'i^.k'i', (1-5) 

where a is the fine structure constant, Cc {rric) is the charge (mass) of the charm quark, 
and k is the energy of the transition photon. The factor I incorporates the matrix 
element for the spin-flip transition of the cc pair. The energy of the transition photon 
is determined by the mass difference between the S-wave spin-triplet and spin-singlet 
charmonium states and therefore is related to the hyperfine mass splitting between 
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them. 

The hyperfine mass sphtting is defined as 

AMMnS) = ^1^(0)^, (1.6) 

where n denotes the principal quantum number, as is the strong couphng constant, 
and V'(O) is the wave function of the S-wave spin-triplet state at the origin. Table [L3] 
lists the theoretical predictions for the hyperfine splitting. The splitting of the first 
radially excited state is consistently less than that for the ground state, reflecting 
coupled-channel effects between the ip{2S) and the nearby open-charm threshold [T^ 
[151 [13]. 

Table 11.41 lists the theoretical predictions for the ip{2S) 7 ric{2S) branching 
fraction. These include phenomenological models, including nonrelativistic potential 
models, effective field theory calculations, and QCD calculations implemented via 
lattice gauge theory (LQCD). The results fall in the range (0.6 — 36.0) x 10^"^. 

1.3.4 Estimates of Production and Decay Rates 

The PDG does not provide enough information for the branching fractions of ip{2S) — > 
^rjc{2S) and ric{2S) —>■ X for definite predictions of event rates in any experiment. 
We have made certain assumptions to estimate these values for our analysis. 

An estimate of expected ip{2S) — > '-yric{2S) rate can be based on the similarity 
between ipC^S) — > 'yric{2S) and J/ip{lS) — > 7?7c(15'). We start from the transition 
rate between spin-0 and spin-1 S-wave states, which is given in Equation II. 5[ For 
estimate, we assume the matrix element is the same for both decays and use the mass 

19 



Table 1.3: Predictions for the hyperfine mass splitting between the first two S- 
wave spin-triplet and spin-singlet charmonium states. The labels AMhf(2S') and 
AMhf(15) are defined as M{i;{2S)) - M(r/c(25)) and M(J/V^) - M{i]c{lS)), re- 
spectively. 



Ref. AMupj^S) [MeV] AMuf{1S) [MeV] 



Experiment [7] 


48 ±4 


116.5±1.2 


Potential Model [16] 


92 


119 


Potential Model [l7j 


83 


130 


Potential Model [18] 


42-65 


99 


Potential Model [19j 


49 


78 


Potential Model [20] 


70 


113 


Potential Model [2 Ij 


72-109 


121-196 


Potential Model [22j 


60 


130 


Perturbative QCD[23J 


69.0 


101.0 


Potential Model 1 241 


71 


116 


Potential Model [25] 


204 


139 


Potential Model [26j 


86 


117 


Potential Model [27j 


78 


117 


HQET[28] 


60 


100 


Potential Model [21] 


68.0 


117.8 


Potential Model [50] 


54.2-82.5 


106.1-128.3 


Potential Model [31] 


98 


117 


BSLT[32] 


38 


102 


Quark Mo del [33] 


98 


117 


Phenomenology [31] 


57±8 




Perturbative QCD[3S] 


37±35 


88±26 


Potential Model [H] 


46.1 


117* 


Quark Model[36] 


51 


119 


Potential Model [37] 


42-53 


108-123 


Potential Model [38j 


38-41 


80-108 


Potential Model [SH] 


66.9-88.4 


115.22-117.06 


Quenched LQCD[l0] 


34±25 


72.6±0.9 


Unquenched LQCD[4l 


-18±189,124±24 


80±1,105±19 



*Input to theory, not a prediction. 
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Table 1.4: Predictions for ip{2S) 7 r]c{2S) branching fraction. The energy of the 
transition photon is denoted by k. The ip{2S) full width (r(ip{2S)) = 337 keV from 
Ref. [7] is used to determine B{ip{2S) — > 7 rjc{2S))^ as described in the text. 



Ref. k [MeV] i3(^(25) ^ 7 r/c(2S)) [IQ-^] 



Experiment [7] 


48±4 




Potential Model [13j 


24.9(49.6) 


0.585(4.45) 


Potential Model [12] 


48* 


3.1 


Potential Model [H] 


92 


0.6-27.0 


Potential Model [22] 


60 


8.0 


Potential Model [13] 


92 


13.4-20.2 


Potential Model [B] 


32 


1.2 


Quark Model[33] 


32(91) 


1.3(29) 


BSLT[15] 


46 


0.89 


Potential Model [37] 


48 


5.0-6.2 


Potential Model [39] 


66.9(88.4) 


12(36) 


Phenomenology [46 ' 


47.8 


2.6±0.7 



* Calculated using mass values from Ref. [7]. 
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differences to calculate the energies of transition photons. For the 25" decay, we have 

E,{i^{2S) ^ 7^c(25)) = ^^'f,. ''^''^ = (47.8 ± 3.9) MeV, (1.7) 

where M^^2S) (-^r,c(25)) is the mass of the ip{2S) = 3686.093 ± 0.034 MeV {ric{2S) = 
3638 ± 4 MeV). Similarly, for the IS" decay, we have 

E^(J/jjj 7770(1-5)) = (114.3 ± 1.1) MeV. (1.8) 

Prom the branching fraction for J/ijj ^ ^rjc, 

BiJ/ij 777e) = (1.3 ± 0.4)% (1.9) 

and the full width of J/-^, 

Tj/^ = 93.4 ±2.1 keV, (1.10) 
we can obtain the partial width of J/-^ — > 777c: 

r(J/V' ^ 7r7c) = B{J/ilj -frjc) x r( J/V') = (1.21 ± 0.37) keV. (1.11) 

With the above information the partial width of ip{2S) jr]c{2S) can be calculated 
as 



/ A-*, I ml ■/>« 1 i. ^/Tl I 1 1 \ ^ 

r(V'(25) ^ 777,(25)) = r{j/i;^^r),{is)) 



= (0.086 ± 0.035) keV (1.12) 

The full width of ilj{2S) is 

rv,(25) = 283 keV. (1.13) 
Therefore, the branching fraction of ip{2S) — > ■jr]c{2S) is estimated to be 

B{ij{2S) 7?7c(25)) = (2.6 ± 1.0) x 10"^ (1.14) 
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To estimate the expected number of produced events for each decay mode con- 
sidered, we need estimates of the branching fractions of r]c{2S) — > X, where X is a 
particular final state. These are estimated by scaling from the corresponding ric{lS) 
branching fractions: 



1.3.5 Previous Measurements of ric{2S) Productions and Decays 

In this analysis we are motivated to measure the mass of the ric{2S) from the decay 
channel 'ip{2S) — 777^(26'), and improve the measurement of full width of the t]c{2S). 
Previous measurements compiled by the PDG give an average of 3638 ± 5 MeV on the 
ric{2S) mass and 14±7 MeV on the full width of ric{2S) in 2005 j^Tj. The measurements 
that contributes to the average of PDG 2006 are from CLEO, M^^(25) = 3642.9 ± 3. 1± 
1.5 MeV and r^^(2S) = 6.3 ± 12.4 ± 4.0 MeV with 61 events in 2004 gZl SSI SH], from 
BaBar, M^^(25) = 3630.8 ± 3.4 ± 1.0 MeV and r^^(2S) = 17.0 ± 8.3 ± 2.5 MeV with 
112 ± 24 events in 2004 [50], and from Belle, M^^(2S) = 3654 ± 6 ± 8 MeV with 
39 ± 11 events in 2002 [51j. In the CLEO measurement, two-photon production of 
r]c{2S), followed by r]c{2S) KsKn was studied with samples 13.6 fb~^ of CLEO 
II/II.V data and 13.1 fb~^ of CLEO III data. The resuh was a measurement of 



the ratio R = [T^^{2S) x B{t]^{2S) KsK7t)]/[T^^{1S) x B{r],{lS) KsKn)] = 



0.18±0.05(stat) ±0.02(syst). The most recent average of measurements gives a small 




(1.8 ±0.9) X B{t]c{1S) X). 



(1.15) 
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correction of the mass of ric{2S) from 3638±5 MeV to 3637±4 MeV, but the width has 
remained unchanged since 2006 [52j. In the BaBar measurement, the r]c{2S) mass and 
width were measured with a sample of 88 fb~^ collected at the T(45'). In 2006, BaBar 
inclusively studied the decay 5^ — > K^ri(.{2S) with a sample of 210.5 fb^^ collected 
at the T(45), finding B{B^ K^t]^{2S)) = (3.4± 1.8(stat) ±0.3(syst)) x 10-^ with 
a significance of 1.8 a [53]. This result is listed in PDG 2008 [52]. In a recent BaBar 
measurement B"^ K^r]c{2S) and t]c{2S) KKtt, where KKn = KsKtt + KKti^ 
was studied with a sample of 349 fb~^ collected at the T{AS). They measured the 
ratio of K^rjc{2S) and ric{2S) KK-n to the corresponding rjcilS) process. 

Using the result from the PDG 2008, this leads to B{t]^{2S) KR-k) = (1.9 ± 
0.4(stat) ± 0.5(syst) ± l.O(br)). However, the mass resolutions of the two separate 
modes a{KsKn) = 9± 1 MeV and a{KKTx^) = 20±9 MeV are very different and the 
peaks of ip{2S) and t]c{2S) overlap each other [54] . In a recent Belle measurement, 
two-photon production of rjc{2S), followed by rjc{2S) Atx or rjc(2S) KKmi, 
was studied with a sample of 395 fb~^ collected at the T(4S'), leading to results of 
T^^{2S) X B{r],{2S) An) < 6.5 eV and T^^{2S) x B{r],{2S) ^ KKtitt) < 5.0 eV at 
90% confidence level [55j . 



1.4 Objectives and Organization of this Dissertation 

Information about ric{2S) production and decay is currently very limited. Measure- 
ments at the T(4S') energy are hampered by serious backgrounds and limited to only 
a few modes. Measurements in the cleaner environment of e'^e~ annihilations near 
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charm threshold have significant advantages. This is the principal motivation for the 
work described in this dissertation. 

Starting in fall 2003, the CESR-c/CLEO-c program began an anticipated three 
year run though the actually run period was 4.5 years until the shut down on March 
2, 2008. During this period, the CLEO-c detector was used to collect the following 
data: 

• 572 pb~^ of integrated luminosity on the ?/'(3770) (see Section [2.1.4p 

• about 27 million ip{2S) decays 

• 21 pb~^ of continuum below the ip{2S) 

• 47 pb~^ of scan data near Ecu = 4170 MeV 

• 13 pb"^ of data at Ecm = 4260 MeV 

• about 600 pb^^ of data at Ecm = 4170 MeV for Ds physics 

While the data samples were smaller than the original project goals, the CLEO-c 
project is regarded as a great success. 

In this project, the decay channel V'(2S') ■yric{2S) is investigated through the 
exclusive reconstruction of candidate events in the CLEO-c ip{2S) data sample. A 
similar nearby decay channel, ^'(25') lXc2^ has been studied for validation and 
comparison purpose, using CLEO's previous measurement of this channel. 

Based on the current information about the decay ipi^S) 'yric{2S), we evalu- 
ated the project by making an estimate on the number of evens that can possibly 
produced. 

25 



For the CLEO-c V'(2S') sample of 25.9 million events, the estimated number of 
ip{2S) 'yric{2S) events produced is A^prod = 6700 ± 2600, where we have used 
a branching fraction of (2.6 ± 1.0) x 10^^. With this information, the expected 
numbers of produced events of all considered modes are listed in Table II. 5[ The 
error in B{ric{2S) —>■ X) is dominated by r{ric{2S))[50%\, and the error in B{iIj{2S) — > 
-fVc{2S)) is dominated by B{J/ij 77/^(25)) [31%]. 

The expected number of produced ric{2S) was estimated assuming a full width 
of 14 ± 7 MeV for r]c{2S) according to PDG average We expected that the 

measurement can be improved with a larger CLEO-c sample based on these PDG 
values as discussed in Chapter [31 

From the above estimate we find that the measurement of the decay ip{2S) 
'jric{2S) is possible in spite of technical challenges and the lack of previous knowledge 
on the physical properties of 'r]c{2S). If enough events can be observed we also plan 
to measure the partial width of 4'{2S) 7?7c(2S'). The only previous contribution 
to the measurement was given by CLEO, in 2004, which is < 0.2 x 10"^ |17l SSJ SS] . 

Following this introduction, the remainder of this dissertation is organized as 
follows. The CESR-c accelerator and CLEO-c detector are described in detail in 
Chapter [21 The analysis set up and the event-selection criteria, as well as the study 
of the related decay channel 4'{2S) — » 7Xc2 are included in Chapter [31 In Chapter [H 
the summary, discussions and conclusions of the investigation of the radiative decay 
4'{2S) 7 rjc{2S) are presented. 
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Table 1.5: Expected number of produced events in each mode. As a reference the 
1% assumption of branching fraction of a mode would provide a number of produced 
event of (67 ±18). 



Mode 


B{ri,{2S) ^ X) 
(%) 


B{^{2S) ^m,{2S))^ 
S(r^,(25)^X)(10-6) 


produced 


47r 


2.2 ± 1.3 


5.7 ±3.3 ±2.2 


147 ± 84 ± 57 


Gtt 


3.6 ±2.3 


9.5 ±5.9 ±3.7 


245 ± 153 ± 94 


KKnn 


2.7 ± 1.8 


7.1 ±4.6 ±2.7 


184 ± 121 ±71 


KKtt^ 


2.1 ± 1.1 


5.5 ±2.9 ±2.1 


142 ± 77 ± 55 


KsKn 


4.2 ±2.2 


10.9 ±5.9 ±4.2 


283 ± 152 ± 109 


7r7r?7(77) 


2.3 ± 1.5 


6.1 ±3.9 ±2.3 


157 ± 100 ±60 


7r7r?7(7r7r7r°) 


1.4 ±0.9 


3.5 ± 2.2 ± 1.3 


91 ± 58 ± 35 


Tmrj', Tj' 7r7rr7(77) 


0.9 ±0.6 


2.3 ±1.5 ±0.9 


59 ± 39 ± 23 


KKrj{'yj) 


< 1.1 


< 2.9 


< 75 




< 0.64 


< 1.7 


< 43 










KKAn 


1.8 ± 1.2 


4.7 ± 3.0 ± 1.9 


121 ±78 ±47 


KsK3n 
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CHAPTER 2 
Experimental Apparatus 



The analysis presented in this dissertation was based on data collected with the 
accelerator and detector at the Cornell University Laboratory of Elementary Particle 
Physics. The design, construction and operation of the apparatus has been carried 
out by the collaboration of hundreds of people. In this chapter the Cornell Electron 
Storage Ring (CESR-c) is described in Section [271] and the CLEO-c detector in Section 
12.21 The trigger and data acquisition systems are described in Section 12.2.71 and 
event reconstruction in Section 12.31 Finally, Section 12.41 explains the simulation of 
the detector response for determination of the detection efficiencies and backgrounds 
for measurements of various physical processes. 

2.1 Cornell Electron Storage Ring - CESR-c 

The Cornell Electron Storage Ring running at charm threshold (CESR-c) is a sym- 
metric electron-positron collider with a circumference of 768 meters. It is located 
at Wilson Synchrotron Laboratory, about forty feet under the track-and-field facility 
of Cornell University. Accelerators like CESR-c produce short-lived particles under 
controlled conditions, allowing more detailed studies than with particles from nat- 
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ural sources like cosmic rays. In CESR-c, electrons and positrons annihilate via 
the electroweak interaction to produce final states consisting of hadrons and leptons. 
These annihilations occur inside the CLEO-c detector, which measures the particles' 
momenta, energies and other properties. 

CESR-c is designed to produce collisions between electrons and positrons with 
center-of-mass energies between 3 and 5 GeV. For the data sample used for this dis- 
sertation, the beam energy was set at around 3.7 GeV, to produce ip{2S) charmonium 
mesons nearly at rest. CESR-c consists of three parts: the linear accelerator (Linac), 
the CESR synchrotron, and the storage ring, as shown in Figure 12.11 

2.1.1 Linac 

The Linac is the first part of the CESR-c, in which electrons and positrons are gener- 
ated, collected, bunched and accelerated before they enter the synchrotron for further 
acceleration. 

The electrons are boiled off of a hot filament by thermionic emission. They are 
accelerated to 150 keV by static electric fields and injected into the 30-meter-long 
vacuum pipe which is the main Linac structure. The Linac has eight radio-frequency 
(RF) cavities, which generate oscillating electric fields that accelerate the electrons 
to an energy of 300 MeV before they are injected into the synchrotron. 

The positrons are created by intercepting the 150 MeV accelerated electrons 
halfway down the Linac and directing them into a tungsten target. Showers of 
electrons, positrons and photons are generated from the target. Positrons are se- 
lected by a magnetic field and then accelerated along the rest of the Linac to a final 
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Figure 2.1: The CESR e+e" collider. 
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energy of 200 MeV for injection into the synchrotron. 

2.1.2 Synchrotron 

The synchrotron is a circular vacuum pipe that fills the CESR-c tunnel and consists 
of four RF accelerating cavities of 3- meter length and 192 bending and focusing 
magnets. The RF cavities add energy to the electrons and positrons each time they 
pass through. Dipole magnets bend the trajectories of the electrons and positrons 
so that they can move along the circular path within the synchrotron. Quadrupole 
and other focusing magnets keep the electrons and positrons confined to trajectories 
near the axis of the beam pipe. As the radius of the circular motion of electrons 
and positrons is given hy R = p/qB, in order to keep the particle beam bunches 
moving in synchrotron, the dipole magnetic field B needs to be synchronized with 
the momentum increase of the particles. It takes about 1500 revolutions or 1/100 
second to accelerate the particles to the desired beam energy, after which they are 
injected into the CESR-c storage ring. At the typical beam energy of 2 GeV , the 
beam particles travel at 99.999997% of the speed of light. 

2.1.3 CESR-c Storage Ring 

Once the electron or positron beam has been accelerated to the desired energy, it is 
injected into the storage ring which operates by the same principle as the synchrotron. 
The electrons and positrons travel in a closed orbit in a much better vacuum than 
the synchrotron. The storage ring contains 106 quadrupole focusing magnets and 86 
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dipole bending magnets, and in addition it also has sextupole and octupole magnets 
for very precise focusing of the beams. When the electrons and positrons circulate in 
the ring, they emit electromagnetic radiation that causes energy loss, which is called 
synchrotron radiation. The energy is restored by superconducting radio-frequency 
(RF) cavities that operate at a frequency of 500 MHz. These RF cavities are similar 
to those used in the Linac and synchrotron, except that they are used primarily to 
restore the energy of the particles, while those in the Linac and synchrotron are used 
to accelerate the particles. 

The synchrotron radiation deposits energy in the vacuum chamber wall and the 
generated heat due to the radiation is carried away by circulating water. Some of 
the synchrotron radiation is used by the Cornell High Energy Synchrotron Source 
(CHESS) facility for X-ray research in the areas of physics, chemistry, biology, envi- 
ronmental science, and materials science. 

In standard operation CESR-c stores the particles in 45 bunches of electrons and 
45 bunches of positrons (configured as nine "trains" of five bunches), circulating in 
opposite directions. They are allowed to cross at one interaction region, where 
the CLEO-c detector is located. One "fill" of electrons and positrons lasts about 
one hour. Since both the electron and positron bunches are stored in the same 
ring, four electrostatic separators are used to set up "pretzel" orbits and ensure that 
the electrons and positrons miss each other when they pass through the unwanted 
intersecting locations, sometimes referred to parasitic crossings. Thus the electrons 
and positrons collide only at the interaction point, which is at the center of CLEO-c 
detector. Figure 12.21 shows a greatly exaggerated [56] schematic of the pretzel orbits. 
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Figure 2.2: Schematic of CESR showing the "pretzel" orbits which are used to 
separate the electron and positron beams at parasitic crossing locations. 
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At the interaction point, the two beams are steered into each other and coUide at a 
small crossing angle of approximately ±2.6 mrad (~ 0.15°). 



2.1.4 Luminosity 

The annihilation rate depends on the "luminosity" of the storage ring, which in turn 
depends primarily on the stored current and beam size. Instantaneous luminosity L 
is defined as the number of collisions for each cm^ of cross section per second by the 
approximate equation 

l^ = nj , (2.1) 

where / is the frequency of revolution of the particles, n is the number of electron or 
positron bunches in the beam, N^- and N^,-^ are the number of electrons and positrons 
in each bunch respectively, and A is the cross-sectional area of the beams. The total 
number of events for a particular process, A^, is given by 

N = ai j Cdt, (2.2) 

where CTj is the cross section for the process. The integral of the instantaneous 
luminosity, J Cdt, is generally referred to as the integrated luminosity or luminosity. 
In CLEO-c, the integrated luminosity for a given data sample is determined by a 
weighted average of the results from three processes e~^e~ — > e~^e~ , and 77, for 

which the cross sections are precisely determined by QED. Each of the three final 
states relies on different components of the detector with different systematic effects 
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2.1.5 CESR-c Upgrades 

CLEO-c and CESR-c have been developed from the previous CLEO III detector and 
CESR accelerator. For over 20 years beginning in 1979, CESR collected electron- 
positron annihilation data at a center-of-mass energy of ~ 10.5 GeV for the study of B 
mesons near threshold. For the CLEO-c and CESR-c project, CESR was proposed 
to run around charm threshold, ~ 3-5 GeV. Modifications to the accelerator and 
detector were necessary to accommodate the lower center-of-mass energy. 

As the rate of the synchrotron radiation is proportional to E"^, the reduction of the 
beam energy greatly altered the beam dynamics. The decrease of the synchrotron ra- 
diation affects CESR performance through two important beam parameters: damping 
time, with which perturbations in beam orbits caused by injection and other transi- 
tions decay away, and horizontal beam size, or horizontal emittance, which measures 
the spread of particles in the bend plane. At decreased center-of-mass energy, these 
effects can limit the luminosity to unacceptable levels. In CESR-c, superconduct- 
ing wiggler magnets were used to increase the synchrotron radiation [56], reduce the 
damping time, and increase the luminosity. A wiggler magnet is a series of dipole 
magnets with high magnetic fields. Each successive dipole has its direction of mag- 
netic field flipped. When a particle passes the wiggler, it will oscillate and therefore 
emit additional synchrotron radiation while keeping the overall path of the particle 
in the ring the same. The damping time is decreased and the horizontal beam size is 
increased because of increase of the synchrotron radiation. CESR-c included twelve 
superconducting wiggler magnets for low-energy running. 
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2.2 CLEO-c Detector 



When an electron and positron collide and annihilate, the flash of energy results in the 
creation of new matter. The CLEO detector, named for the one historically coupled 
to Caesar, was built to study these collisions, with details described in References [58| 
[59]. 

The CLEO-c detector was the final stage of a series of upgrades of the CLEO 
detector since its commissioning in 1979. Beginning in 1989, the CLEO II detector 
was a productive source of physics research, due to its very good drift chamber based 
tracking system and outstanding electromagnetic calorimeter. In 1995, after a three- 
layer silicon strip vertex detector and a new beam-pipe were installed, and an argon- 
ethane gas mixture in the drift chambers was switched to helium-propane, CLEO II 
was upgraded to CLEO II. V. A few years later CLEO II. V was upgraded to CLEO 
III by adding a Ring Image Cherenkov (RICH) detector for particle identification and 
a new tracking system, which consisted of a new central drift chamber surrounding a 
four-layer silicon strip vertex detector. The conversion from CLEO III to CLEO-c in 
2003 consisted of the replacement of the silicon vertex detector with a low-material 
gaseous tracking detector (ZD). 

2.2.1 CLEO-c Detector 

The CLEO-c detector (Figure 12. 3p measures about 6 meters on a side, and is con- 
structed of about 900,000 kilograms of iron and over 25,000 individual detection 
elements. It is an approximately cylindrically symmetric detector aligned along the 
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beam line, and covering about 93% of the 47r solid angle. Electrons and positrons 
collide at the interaction point (IP) in the beam pipe, which is located at the center 
of the detector. The annihilation produces new particles, which decay quickly into 
long-lived or stable charged particles, electrons, muons, protons, pions, and kaons, 
and some neutral particles including photons. Particles pass through different sub- 
detectors of the CLEO-c detector, leaving measurable signals. The tracking system 
provides the precise measurement of the trajectories of charged particles and infor- 
mation on the rate of energy deposit in the material (dE/dx). The Ring Imaging 
Cherenkov (RICH) detector provides a velocity measurement for charged particle 
identification. The electromagnetic crystal calorimeter (CC) is used to measure the 
energies of electrons and photons. A superconducting solenoid provides a uniform 
magnetic field over the whole tracking system, RICH and CC. The muon detector 
is used to detect muons, which can penetrate all inner sub-detectors if they have 
sufficient momentum. 

The whole operation of the experiment includes the constant calibrations, event 
data taking process and offiine data analysis, which includes the Monte Carlo studies 
and physics analysis. 

2.2.2 Tracking System 

The tracking system of CLEO-c consists of two drift chambers, which provide precisely 
measured space points for charged particles emerging from the interaction. These 
devices reside inside a superconducting solenoidal magnet that provides a uniform 
field of 1 T along the beam direction. The tracking system measures the helical 
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Figure 2.3: The CLEO-c detector. 
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trajectories of charged particles, and online and offline software reconstruct these 
trajectories to obtain precise measurements of particle momenta. 

Inner Drift Chamber - ZD 

The innermost sub-detector of the CLEO-c detector is the inner drift chamber (ZD), 
which is located right outside of the beam pipe. The major hardware modification 
from CLEO 111 to CLEO-c is the replacement of the Silicon Vertex Detector with 
the wire vertex drift chamber ZD. The silicon vertex chamber, as the tracking drift 
chamber for CLEO 111, was used to provide extremely accurate track position mea- 
surements in r — {/) and z as close as possible to the interaction point for studies at the 
T(4S'). However, when the center of mass energy was reduced from 10 GeV down to 
3-5 GeV, the momentum distribution of the charged tracks is shifted down to lower 
values. Minimizing material is crucial at lower center-of-mass energies since multiple 
scattering dominates the momentum measurement error for low momentum tracks. 
As the tracks have lower velocities, the importance of vertexing is reduced because 
vertex separation for decays of particles like and D° are too small to be resolved 
with the silicon vertex detector. Since the extreme precision of the previous silicon 
vertex detector was not needed and the material was a disadvantage, it was replaced 
with the ZD chamber [50] . 

The ZD is located just outside the beam pipe and covers radii from 4.1 cm to 
11.8 cm from the interaction point. An isometric view of the structure of the ZD is 
shown in Figure 12.41 It supports gold-plated tungsten sense wires and gold-plated 
aluminum field wires. The sense wires are in six layers with the number of cells per 
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Figure 2.4: An schematic of ZD drift chamber. 
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layer ranging from 34 to 66 for a total of 300. Each cell measures 10 mm and consists 
of a sense wire held at +1900 V, surrounded by field wires at ground. Cells share 
bordering field wires. The wires are oriented at a small (stereo) angle with respect 
to the beam axis, allowing measurement of r, and z. The stereo angle of the ZD 
ranges from 10.3° in the innermost field wire layer to 15.4° in the outermost field wire 
layer, with the stereo hyperbolic sag of 4 mm. The "all stereo" design optimizes 
the measurement of z (position in the beam direction), from which the inner drift 
chamber was given the name "Z Detector." The ZD's stepped structure provides 
spatial measurement of charged particles within | cos 6*1 < 0.93, where Q is the angle 
of the the particle with respect to the beam. 

When a charged particle passes through a cell, it ionizes the atoms of gas filling. 
The ZD is filled with a very pure mixture of 60% helium and 40% propane (CsHg) 
which provides a very long radiation length (~ 330 m) [56j. The electrons produced 
from the ionized atoms are accelerated toward the sense wire and away from the field 
wires. As the electric field increases toward the sense wire, these primary electrons 
ionize other gas atoms and create an "avalanche" at the sense wire. The time of 
the electric pulse observed on the sense wire and the charge, which is related to the 
deposited energy, are recorded. This information is used to map out the trajectories 
of the charged particles and to fit tracks. 

Main Drift Chamber - DR 

Outside of the ZD detector is the main drift chamber (DR), a much larger device 
that measures particles out to a distance of 0.8 m from the interaction point. The 
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main drift chamber has design similar to that of the ZD, but with much larger size. 
It has 9,796 gold-plated tungsten sense wires and 29,682 gold-plated aluminum field 
wires, arranged in 47 layers. The wires are grouped in nearly square cells, with cell 
size 14 mm. Each cell contains one sense wire surrounded by eight field wires. The 
voltage on the sense wires is -1-2100 V with respect to the grounded field wires. The 
gas in the DR is the same as the gas in the ZD, a mixture of 60% helium and 40% 
propane. Among the 47 layers of wires, the first 16 are axial and the remainder are 
stereo layers that alternate angles about 3° in groups of four to provide z information 
throughout the volume of the detector. In order to determine the z position of 
the particle at the stereo layer, the axial layers are used to predict the particle's 
r — (f) position and to match the r — (f) information from the stereo wires. The inner 
and outer radii of the drift chamber are covered with longitudinally and azimuthally 
segmented cathodes to provide precise measurements of the z position of most tracks 
as they enter and emerge from the DR. 

Hadron Identification by Specific Ionization - dE/dx 

Information about the rate of energy loss is provided by the DR for each charged 
particle, along with the the measured momentum. This can be used for particle 
identification. The energy lost per unit length depends upon a particle's velocity, 
as given by the Bethe-Bloch formula. A variable is formed for each particle 
hypothesis of electron, muon, pion, kaon, or proton. The value of is the sum of 
Xi over a lot of hits. Xi is defined for hit i as 

_ ciE'/da; (measured) — dE / dxi{expected) 
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where Ui is the uncertainty of the measurement, usually approximately 6%. Figure [231 
shows the measured dE/dx as a function of momentum for a large population of 
charged particles. At momenta below 500 MeV/c, pions and kaons are well separated. 
At higher momenta, the dE/dx information is limited and additional information is 
needed for good particle identification. 

Track Reconstruction 

Track reconstruction is the process to form tracks with hits, the electrical signals from 
the wires, using a pattern-recognition algorithm on three-dimensional (r — and z) 
position measurement from the ZD and DR. (The r — cj) resolution varies over the 
cell and is of order 100 /im.) Once tracks are found they can be fit to determine 
interesting physics quantities like momenta, vertex positions and directions [61]. The 
CLEO-c fitter is an implementation of the Billoir or Kalman algorithm that optimizes 
the determination of the momentum and trajectory, taking into account the expected 
energy loss. 

The tracking system is within an axial magnetic field of 1 T. The charged particles 
follow helical paths within the constant magnetic field. As the magnetic field is 
parallel to the z axis, the transverse momentum Pi, in the r — cj) plane, can be related 
to the curvature of a particle's trajectory by 

= qBa, (2.4) 

where q is the magnitude of the particle's charge, B is the magnitude of the magnetic 
field, and a is the radius of curvature. Therefore, the measurement of the full 
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Figure 2.5: Measured dE/dx in the CLEO-c DR as a function of measured momen- 
tum. The bands represent particles of different mass: electrons, pions and kaons 
from left to right. 
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momentum vector can be obtained by the measurement of the track's curvature in 
r0, and its angular coordinates 9 and 0. At 1 GeV/ c the momentum resolution of the 
charged particle is about 0.6%. The direction of the curvature within the magnetic 
field indicates the sign of the charged particle. 

2.2.3 RICH Detector 

The Ring Image Cherenkov Detector (RICH) [621 [63] is located directly outside of 
the main drift chamber. It was the most significant improvement in the CLEO III 
upgrade from CLEO II. V and was not modified for CLEO-c. The measurement of 
Cherenkov radiation in the RICH provides the additional charged particle identifica- 
tion needed for the higher momenta where dE/dx is inadequate. 

When a charged particle passes through a transparent dielectric medium at a 
speed greater than the speed of light in that medium, electromagnetic radiation, 
called Cherenkov radiation is produced. This effect was first observed by Cherenkov 
in 1934. As a charged particle travels through the medium, its electromagnetic field 
disturbs the local electromagnetic field in the medium, and displaces and polarizes 
the electrons in the atoms of the medium. Photons are emitted when the electrons 
of the medium return to their equilibrium state. When the charged particle travels 
at a speed higher than the speed of light in the medium, the photons constructively 
interfere to make the radiation observable. The radiation is produced in a cone with 
its central axis along the trajectory of the traveling particle. The characteristic angle 
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of the cone 6c, known as Cherenkov angle, is related to the velocity of the particle by 



where n is the index of refraction of the dielectric medium, and f3 is the ratio between 
the velocity of the charged particle and the speed of light c. With the relations 
(3 = p/E and E"^ = m? + the cosine of the Cherenkov angle can be expressed in 
terms of index of diffraction n, mass of the charged particle m, and momentum of the 
charged particle p by 



If the momentum of the charged particle can be measured independently, such as by 
fitting the trajectory in the tracking system, then the Cherenkov angle gives the mass 
and therefore the identity of the charged particle. 

The RICH detector is shown in Figure 12.61 and covers 83% of the 47r solid angle. 
The innermost parts of the RICH detector are the radiators, which are made from LiF 
and have an average thickness of 1.7 cm. For CLEO-c's magnetic field of 1 T, tracks 
with transverse momenta of 0.12 GeV/c reach the RICH radiators. When a charged 
track passes through the radiators, Cherenkov photons are produced. There are a 
total of 14 rows of radiators and the four central rows contain "sawtooth" radiators, 
with triangular grooves on the surface, to overcome the total internal reflection for 
Cherenkov photons produced with nearly normal track incidence. The outer rows 
of radiators are fiat. The Cherenkov photons produced in the radiators propagate 
outward through an expansion volume filled with pure nitrogen, which is transparent 
to Cherenkov photons with the characteristic wavelength of about 150 nm. Prop- 




(2.5) 
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Figure 2.6: The r — <j) cross-section view of the CLEO-c RICH detector. 
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agation through the expansion volume allows the cone of radiation to be detected 
where its radius is large enough for accurate measurements. After expansion, the 
photons pass through CaF2 windows and enter the multi-wire proportional chambers 
(MWPC). The chambers are filled with a gas mixture of methane (CH4) and tri- 
ethylamine (TEA, N(CH2CH3)3) [60], which converts Cherenkov photons in a narrow 
ultraviolet range of 135-165 nm to photo-electrons. Gas multiplication of the pri- 
mary ionization occurs as described for the ZD and DR, although in the RICH the 
location of the photon conversion is measured through charge induced on cathode 
pads. There are 230,400 pads in total. Sample Cherenkov ring images are shown in 
Figure O 




Figure 2.7: Cherenkov ring images produced by a charged track crossing sawtooth 
(left) and flat (right) radiators. The rectangular grids are partial 24 x 40 pad arrays 
and the small squares represent the charge detected in pads. The hits at the center 
of the ring are due to the charged track crossing the wire chamber. The other 
hits are due to produced Cherenkov photons. The image of the flat radiator only 
shows half of the Cherenkov ring as the other half is trapped in the radiator by total 
internal reflection. The sawtooth radiator image is distorted by refraction. 
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To identify particles the radius of the Cherenkov cone is measured and combined 
with the expectations for a particular hypothesis to construct a likelihood ratio. For 
two different particle hypotheses, a difference variable for likelihood ratio is rep- 
resented by 

X?-x| = -21nL, + 21nL,, (2.7) 

where i, j are different particle hypotheses and can be electrons, muons, pions, kaons, 
or protons, and Lj, Lj are their corresponding likelihoods. Figure 12.81 shows the 
measured fraction of kaons and pions as a function of the cut on — x^. When the 
cut is set at xl: ~ < 0; 92% of kaons can be identified with a pion fake rate of 8%. 
The measured momentum range of kaons or pions is from 700 MeV/c, which is just 
above kaon radiation threshold, to 2.7 GeV/ c. Figure [Z9] shows the separation of the 
different particle hypotheses as a function of momentum. The curves are plotted with 
the momenta of particles for both particle hypotheses to be above their respective 
thresholds, which is determined by the index of refraction of the LiF radiator, n = 1.4. 

For momenta below the threshold, RICH is used in a threshold mode and the 
number of photons assigned to the pion hypothesis is used instead of Xk ~ xi ■ The 
kaon efficiency and pion fake rate as a function of the number of photons assigned to 
the pion hypothesis are shown in Figure 12.101 

2.2.4 Electromagnetic Calorimeter 

The electromagnetic crystal calorimeter (CC) is right outside of the RICH detector 
and inside of the superconducting solenoid, as shown in Figure 12.31 and 12.111 It is 
used to measure the energies of electrons and photons, and covers approximately 93% 
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Figure 2.8: Kaon identification efficiency (filled circle) and pion fake rate (open cir- 
cle) as a function of various cuts on difference between kaon and pion hypotheses. 
The momentum of tracks is between 0.7 GeV/c, just above kaon radiation threshold, 
and 2.7 GeV/c. 
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Figure 2.9: Separation of particles with different particle hypotheses in the RICH 
detector as a function of momentum. All curves have cuts at minimum momentum 
because both particles are required to be above their radiation threshold, determined 
from the index of refraction of the LiF radiator, n = 1.4. 
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Figure 2.10: Kaon identification efficiency (filled circles) and pion fake rate (open 
circles) as a function of the required number of photons assigned to the pion hy- 
pothesis. The momentum of tracks is less than 0.6 GeV/c. The kaon radiation 
threshold is 0.44 GeV/c. 
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Figure 2.11: r — z cross-section view of CLEO-c detector. 
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of the full 4:71 solid angle. It consists of 7784 thallium-doped Cesium Iodide (Csl) 
scintillation crystals. There are 6144 crystals arranged in the barrel region, defined 
by |cos^^| < 0.80, in a projective geometry pointed roughly toward the interaction 
point. The two endcap regions have 1640 crystals, covering 0.85 < |cos^^| < 0.93. 
Each crystal measures approximately 5 cm x 5 cm square by 30 cm long. The energy 
resolution of the barrel and endcap regions are slightly different, as illustrated in 
Figure 12.121 The energy resolution of the transition region, between the barrel and 
endcap (0.85 < |cos^^| < 0.93) is degraded by the substantially larger amount of 
material in front of the CC in this region. It is excluded in most analyses. 

The material from which the crystals were made, Csl, has high density (4.53 g/cm^) 
and short radiation length (8.93 g/cm^), so it can stop photons and electrons effec- 
tively. In the calorimeter, electrons produce photons by Bremsstrahlung and photons 
undergo e+e~ production near high-Z nuclei. This process produces a cascade called 
an electromagnetic shower. The final low-energy electrons and positrons that result 
from this shower process excite atoms in the crystals which then "scintillate," emitting 
visible light (560 nm) as they return to the ground state. The crystals are transparent 
for these visible light photons and they are read out using four 1 cm x 1 cm pho- 
todiodes mounted on the backs of the crystals. Muons and noninteracting hadrons 
are called minimum ionizing particles (MIPS) as they only deposit a small fraction 
of their energy by ionization. 

The shower reconstruction process starts with calculating the energy deposited in 
a crystal from the amount of light detected. As the energy of a shower is generally 
deposited into multiple crystals, a cluster of several adjacent and near- adjacent crys- 
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Figure 2.12: M^j resolution for 7r° candidates in CLEO III T(45) data. 
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tals are considered together. Since including crystals with very small energy deposits 
can degrade the measurement by bringing in excessive noise, the number of blocks 
used has been optimized for the best precision. The uncertainty in the energy of 
the shower is calculated based on energy-weighted average of the blocks in the sum 
of the cluster, and an example optimization is shown in Figure I2.13[ The number of 
crystals used is logarithmic in energy and ranges from 4 at 25 MeV to 17 at 4 GeV 
[58]. 

The photon energy resolution of CLEO-c calorimeter is about 4.0% at 0.1 GeV 
and 2.2% at 1 GeV and the angular resolution is approximately 10 mrad. The per- 
formance of the endcap region is not as good as that of barrel region because of the 
presence of the aluminum DR endplates and electronics in front of the crystals. The 
excellent energy and angular resolution of the calorimeter is critical for the recon- 
struction of the decay modes that include vr'' ^ 77 or 77 — > 77 decays, as well as the 
low-energy transition photons in ip{2S) radiative decays. 

2.2.5 Superconducting Magnet 

Right outside of the calorimeter is the superconducting solenoid, which has an inner 
diameter of 3 m and a length of 3.5 m. It provides a uniform (to ±2%) 1-T mag- 
netic field parallel to the beam line over the full volume of the tracking system, the 
RICH detector and the calorimeter. The solenoid is cooled down to superconducting 
temperatures by liquid helium. Three layers of 36 cm thick iron flux return for the 
magnet also serve as part of the absorber for the muon identiflcation system. The 
magnetic field used to be set at 1.5 T for CLEO III running at a center-of-mass energy 
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Figure 2.13: Energy resolution as a function of the number of crystals (blocks) 
used for shower reconstruction in CLEO II. The same calorimeter has been used 
from CLEO II to CLEO-c. The smooth lines represent Monte Carlo simulations of 
100 McV photons with different noise levels. The points are the data from 100 McV 
transition photons of the decay T(3S') 7x^7 (2-P), and the arrow indicates the 
number of crystals used for the 100 MeV photon reconstruction. 
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of 10.5 GeV. In CLEO-c the center of mass energy is much lower and the momentum 
of charged particles is much smaller. From Equation 12.41 one can see that in order to 
keep the same radius of the curvature for the trajectory of particle, if the momentum 
of the particle decreases, the magnetic field has to be reduced to allow the tracks to 
reach the RICH detector and calorimeter. 

2.2.6 Muon Detector 

Muons are distinguished from other charged particles by their low probability of 
interaction as they pass through material. The Muon Detector (MU) is located at 
the outside of CLEO-c beyond a roughly 1-m steel absorber. The absorber absorbs 
other particles that pass through the calorimeter, so the hits in the MU chambers are 
almost exclusively from muons. The muon chambers (MU) operate in the similar 
way the tracking chambers do. A cross-section view of the muon super-layer is shown 
in Figure [2.141 The muon detector uses gas-filled tracking chambers in between 36 cm 
iron absorbers surrounding the detector. It consists of three layers of proportional 
wire chambers of 4 m long, 8.3 cm wide, and 1.0 cm tall and oriented with its long 
axis parallel to the z-axis. 

The muon detector provides the information in terms of interaction length. The 
interaction lengths (also known as the absorption length) is defined as the mean 
distance over which a particle travels before scattering inelastically from a nucleus. 
The barrel super layers are located at depths of 3, 5, and 7 interaction lengths, while 
the only super-layer in each endcap is located at a depth of about 7 interaction 
lengths. 
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Figure 2.14: A cross-section view of a muon chamber super-layer. 

The MU is designed to detect muons with momentum of 0.8 Gev/c or higher. It 
is not as useful for CLEO-c as it was for the previous experiment because of the much 
lower center of mass energies. 

2.2.7 Triggers and Data Acquisition System 

The trigger and data acquisition (DAQ) systems of CLEO-c are used to collect and 
process the data from all components of the detector, and save the data to mass 
storage devices. During CLEO-c runs, interactions happened at the frequency in 
the order of 1 MHz. The recording capability of the data acquisition system does 
not allow every collision to be recorded, and most of these interactions are of no 
interest. The actual rate of interesting physics events is only on the order of 1 Hz. 
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The trigger system, as illustrated in Figure I2.15[ is used to decide whether an event 
should be recorded based on the complicated information from many components of 
the detector. 

Trigger System 

The trigger system is operated in several levels [561 EH [651 ES]. The first level of 
trigger, called level- zero (LO), is the hardware trigger. The tracking trigger, includ- 
ing axial and stereo tracking, and calorimeter trigger, including analog and digital 
calorimeter trigger, both belong to LO. Data from the drift chamber and calorime- 
ter are received and processed in separate VME crates using the appropriate circuit 
boards to yield initial trigger parameters, such as track count and topology in the drift 
chamber, and shower count and topology in the electromagnetic calorimeter. The 
next level trigger is the global level-one (LI) trigger, which correlates the information 
from both tracking and calorimeter systems. The global trigger circuit generates an 
LlPass strobe every time a valid trigger condition is satisfied. The generated LlPass 
signals are conditionally transferred to the gating and calibration modules and then 
distributed to the data acquisition system. The LI Decision and Data Flow Control 
system makes a trigger decision every 42 ns based on the information from the track- 
ing system and calorimeter. Programmable trigger decision boards (LITR), which 
can be configured for a wide range of topologies, are used to look at this informa- 
tion and respond to various trigger conditions. It takes the tracking system about 
2 fis and calorimetry over 2.5 fis to deliver information for the decision. There are 
currently about eight trigger lines, which definitions and relative rates are shown in 
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Figure 2.15: Overview of the trigger system. 
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Table [XTl The total trigger rate at a luminosity of £ = 1 x 10'^^ cm ^ is between 
80 and 90 Hz [56]. 
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Table 2.1: Definition of CLEO 111 trigger lines. The rate is relative to barrel 
Bhabhas line. 

The final trigger stage is called level-three (L3). It is a software trigger, which is 
implemented on a fast workstation. It takes an event passed by the lower level triggers 
and either rejects or accepts it based on more sophisticated reconstruction algorithms. 
Events such as those from cosmic rays or interactions of the beam with residual gas 
molecules in the beam pipe are rejected by L3 before they are permanently recorded. 
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Data Acquisition System 

The data acquisition system consists of two equally important parts [56] . Data qual- 
ity and the detector performance are monitored by the slow control system, and the 
data collection system is responsible for the data transfer from the front end electron- 
ics to the mass storage device. The slow control and monitoring system is required 
to provide sufficient flexibility for online data quality and detector performance mon- 
itoring. It also includes run control, alarm and message handling, and calibration 
constants management. For the data collection system the knowledge of cross sec- 
tions of the interesting physics processes and expected CESR luminosity are needed 
to establish the performance requirement. The key parameters for the DAQ system 
are the expected trigger rate, the acceptable detector dead-time and the average event 
size. In the current configuration the data acquisition system can operate at trigger 
rates up to 150 Hz, significantly below the design rate of 1000 Hz. A read-out time, 
defined as the time between the trigger signal and the end of the digitization process 
in the front-end electronic modules, causes dead-time and therefore should be kept 
minimal. The detector dead-time is less than 3% when the average read-out time is 
20-30 fis at a trigger rate of 1000 Hz. 

The structure of the data acquisition system is shown in Figure 12.161 For each 
event that is accepted by the trigger, all detector channels need to be digitized. Front- 
end data conversion is performed in parallel and local buffers on each data-board 
hold the data for later asynchronous readout by the data acquisition system. Data 
sparsification (the suppression of channels without valid hits) is performed directly 
on the data-boards. The Data-Mover, a dedicated module in each front-end crate, 
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Figure 2.16: Architecture of the CLEO III Data Acquisition System. The only 
change from CLEO III to CLEO-c in this diagram is to replace the silicon vertex 
detector by the ZD. 
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assures transfer times below 500 and provides a second buffer level. About 30 
front-end crates are needed for the detector. The data is sent to L3 and the L3 
decision is sent back to the front-end crates. The fragments of the events that are 
accepted by L3 are transmitted from the crates to the Event-Builder by the Data- 
Mover. Completely reconstructed events are transferred to the mass storage device 
and a fraction of the data is analyzed online by a monitor program to quickly discover 
problems of the detector and to ensure the quality of the stored data. 

2.3 Event Reconstruction 

Before the data can be used for physics analyses, the raw data collected needs to 
be processed by software packages to reconstruct tracks and showers. Raw data 
collected from different components of the detector are processed through Pas si and 
Pass2, which consist of collections of processors and routines for event reconstruction. 
The detailed processes of the track and shower reconstruction have been described in 
Sections [2:2:21 and [2:2:1 

The function of Pas si is to perform fast classification of the raw data events. 
It runs online to provide real-time feedback about the performance of the detector 
while taking data for data quality control. Only one-tenth of the collected events are 
processed with fast reconstruction. Pass2 processes all events offline and it provides 
more detailed and accurate event reconstruction. After Pass2 the data are ready for 
the physics analyses and one can access the data with personal computer program 
written with a computer language like C+-I-. 
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2.4 Monte Carlo Simulation 



Monte Carlo (MC) techniques allow the simulation of the data based on the best 
knowledge of physics properties, detector configuration and detector performance. 
The simulated data are used to determine detection efficiencies, estimate backgrounds, 
or determine background and signal shapes for specific physical processes studied, 
including in this dissertation. 

The CLEO-c MC production [671 starts with the generation of a list of particles. 
In this stage, for example, 'ip{2S) events are generated with the software package 
EvtGen [67] and final state radiation (FSR) is simulated with PHOTOS [69]. EvtGen 
produces Monte Carlo events with a list of particles (and their four-momenta) created 
in an e+e~ interaction at a desired energy. In addition to specific states like ■iIj{2S), it 
simulates events for e~^e~ qq where q can be an up, down, strange, or charm quark. 
The quark fragmentation simulation uses an interface to the JetSet software package 
[70l 171] that models quark and gluon hadronization. The subsequent decays of 
unstable particles are simulated based on the amplitudes, which are used to calculate 
the probability for any decay tree from a decay table that is updated periodically to 
reflect our best knowledge. 

The CLEOg is the simulation program based on the GEANT3 [72j package de- 
veloped at CERN. CLEOg simulates the passage of particles through the CLEO-c 
detector material and detector elements. At each step, a random number is gener- 
ated to determine if the particle interacts with any material. If it does, more random 
numbers are generated to determine its effect on the particle and still more to cal- 
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culate the signal that would result in the detector element. Using calibration data 
or samples of data events taken with a random trigger, the detector noise and the 
beam-related backgrounds are also simulated when no colliding physics events were 
present. Finally, the simulated events are saved in the same format as raw data so 
they are processed with the same reconstruction program Pass2. 
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CHAPTER 3 

Measurement of '0(25') 777c (25) 



In this chapter, the analysis of the process of ip{2S) 7?7c(25') is described in detaiL 
As the channel had never been accurately measured before, I start by introducing 
the strategy and procedures of this analysis, followed by describing the samples used 
and the event selection criteria. The study of ip{2S) — » 7Xc2 decays is presented to 
justify the analysis procedures of ipi^S) ■yric{2S) decays in Section [3751 Finally the 
results of i'{2S) — > ■yric{2S) decays and systematic uncertainty studies are presented. 



3.1 Strategy and Procedures 

The analysis was initiated from the CLEO exclusive study of 4'{2S) — > jric{lS), 
which was included in the search for ip{2S) vr^/ic [731 [Tlj. We reconstruct 4'{2S) — > 
^r)c{2S) candidates and extract information from the energy of the transition photons. 

Various modes have been searched. We constructed an extensive list of possible 
decay modes following the model of a CLEO-c analysis of J/ip —>■ 'yrjci^S) [75]. We 
considered modes with relatively high r]c{lS) yields, and selected the following final 
states: An, Gvr, KKtttt, KKtt^, KsKtt, -k-kt] with — ^ 77 or -r-kti^ , ttttt]' with rj' 
7r7r?7(77), KKt] with 77 77 or tttttt^, KKinrT:^, KKAtt, and KsKStt. 
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The major challenge of this analysis is to detect the low energy transition photon 
while suppressing backgrounds. We reconstruct the events exclusively by considering 
all combinations of charged tracks and photons, and remove the possible candidates 
of other decays by applying various mode dependent cuts, as described in Section 
13.31 Signal- squared over signal-plus-background [S'^/{S + B)] studies were used to 
optimize selection criteria for the 4-C kinematic fit x^/dof and the distance to the 
nearest track for the transition photon and the angle between the transition photon 
and the closest track. 

Some other modes have been considered but were not included in the final analysis. 
The modes with large ip{2S) decay branching fractions are expected to have large 
backgrounds. These modes include 47r27r° and KsKst^ti- In addition to these 
two modes, the following modes are not considered because of high background level 
based on studies of generic ip{2S) Monte Carlo (MC): nnr]', r( — > 7p°(7r7r), 27r27r°. The 
modes Trvrr/',//' 7r7r?7(7r7r7r°), 47^^, KsKtttt^ are not considered due to poor signal- 
to-background ratio based on small expected yields from extrapolating the partial 
widths from the ?7c(l>S') to the ric{2S). Modes KKAtttt^ and KsKSttti^ are excluded 
because signal MC studies show they have low efficiency. When optimizing the cuts 
for the considered modes, the branching fraction of each ric{2S) mode was arbitrarily 
assumed to be 1% because there are no previous measurements that can provide 
reliable estimates. 
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3.2 Samples 



3.2.1 Data 

In order to search for ip{2S) — > jric{2S), CLEO-c Datasets 32 and 42 were used. 
Dataset 32 has 1.44 milhon 'ip{2S) decays and Dataset 42 has 24.45 milhon, giving a 
total for the CLEO-c data sample of 25.89 million ip{2S) decays [76]. For our final 
processing we selected a recent update of the CLEO-c software that allowed us to 
use updated endcap calorimeter (CC) calibrations and a new CC energy-correction 
procedure (CCFIX) [TT]. 

3.2.2 Generic Monte Carlo 

The 5 times luminosity generic il^{2S) and the 5 times luminosity continuum MC 
samples, both of which were generated by the Minnesota MC farm, were used for 
cut optimization and background suppression. The full 10 times generic ip{2S) and 
continuum MC samples were used for fitting the measured photon energy distribution 
in the ric{2S) signal region. The new 5 times luminosity generic ip{2S) sample is found 
to be consistent with that of the original 5 times luminosity sample for this study. 

3.2.3 Signal Monte Carlo 

The signal MC sample consists of one million events, distributed over 15 modes, 
including the 13 modes that were used in this study. (We count the modes with an 
Tj decay as two separate modes, although they are combined for final measurements.) 



70 



The sample was generated to replicate the full CLEO-c ipi'^S) data sample. The 
r]c{2S) mass (3638 ± 4 MeV) and width (14 ± 7 MeV) were taken from the PDG 
This MC sample was generated using the appropriate angular distribution for a 
vector to vector-pseudoscalar decay, i.e., P(cos6') = 1 + cos^6', where 9 is the angle 
between the Ml transition photon and the positron beam. The angular distributions 
for the ric{2S) decays were thrown according to phase space. The minimum hadronic 
mass of the ric(2S) was decreased from its default value of 20 MeV below the nominal 
mass of ric{2S) (3638 MeV) to 210 MeV below in order to appropriately model the 
tails of the r]c{2S) decay. 

A signal MC sample to study the Xc2 decaying to the same hadronic final states was 
also generated. The study of Xc2 decays was used to validate the analysis procedure. 
A one million event sample, distributed over the same modes, was generated. This 
Xc2 MC sample was generated using the appropriate angular distribution for a vector 
to axial vector-pseudoscalar decay, i.e., P(cos6') = 1 + (1/13) ■ cos^ ^ [78]. As with 
the i]c{2S) signal MC sample, this sample was generated to replicate the full CLEO-c 
ip{2S) data sample. The Xc2 mass (3556.20 MeV) and width (2.05 MeV) [7] were 
used. The minimum hadronic mass of the Xc2 was decreased from its default value of 
6.0 MeV below the nominal mass of Xc2 (3556.20 MeV) to 41.6 MeV below in order 
to appropriately model the tails of the Xc2 decay. 

The numbers of generated events for both 7]c{2S) and Xc2 decays are listed in 
Table O 
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Table 3.1: Number of generated signal MC events. The resonance parameters used 
to generate these samples are described in the text. 



Mode 






47r 
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71388 


Qtt 
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KKtt^ 
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71248 


KsKtt 


73383 


70853 


TTTT 77(77) 


71052 
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7r7r?7(7r7r7r°) 


70928 
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TTTrr]', Tj' — > 777777(77) 


70986 


49970 
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71145 




71261 
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KKnnn'^ 
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KKAn 


62426 


70908 


KsKStt 


64556 


71063 
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3.3 Event Selection 



It was required that the net charge of the event equal zero and the number of good 
tracks match the number of charged tracks in the decay mode. 

3.3.1 Track Selection 

We used standard CLEO-c selection criteria for track quality and particle identifica- 
tion (PID). These have been developed for use in high-precision measurements of D 
decays by the CLEO-c "DTag" subgroup [79] ■ They are applied in selecting K^'s 
and TT^'s. 

Track Quality 

For track quality, it was required that the track have a successful Kalman fit. Also 
the projection of the distance from the interaction point of the e^e~ annihilation (IP) 
to the origin of the track on the r — cj) plane \Dq\ < 5 mm and the projection of the 
distance between the IP and origin of the track along the z-axis \Zq\ < 5 cm were 
required for each qualified track, with exception of charged tracks from Ks tt^tt^ 
decays. (If charged tracks are daughters of a Ks, the Dq and Zq of the track tend 
to be larger as the charged particles are created at where the Ks decays rather than 
the IP.) 
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Particle ID 

The Particle ID criteria developed by the DTag group were used for charged pions 
and kaons. Those criteria are 

• I cos^l < 0.93, where 9 is the angle between the initial momentum of the track 
and the direction of the beam, the z axis 

• RICH information is valid 

• Both TT^, hypotheses analyzed 

• 3(j dE/dx consistency with the vr^ or hypothesis, which means if the dE/dx 
distribution for the vr^ or hypothesis is Gaussian, then 

I ((i£;/(ix)measured - {dE / dx)^± K±\ ^ g 
<^ dE/dx 

where {dE/dx).,^± or k± and (J dE/dx are the mean and standard deviation of the 
the dE/dx distribution for the vr^ or hypothesis 

• If the momentum of the particle is greater than or equal to 700 MeV: 

o More than two photons on the Cherenkov ring (See Section 12. 2. 31) are 
associated to the pion or kaon 

o Combined log-likelihood L = Lt,± — Lk + cr^± — > for tt^, < 
for K^, where L^^i and L^± are the log-likelihoods for the vr^ and K"^ 
hypothesis, respectively, and cr^± and aK± are the standard deviation of 
the dE/dx distribution for the vr^ and hypothesis as defined above as 

(^dE/dx 
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• If the momentum of the particle is smaller than 700 MeV, or if RICH information 

is not available: a^± — cr^± > for tt^, < for 

For any mode including a K^K~ pair, only one kaon is required to pass the PID 
to increase the efficiency. 

3.3.2 Photon Selection 

General requirements for selecting photons, both from radiative transitions and in 
7r°, ?7 ^ 77 reconstruction, were applied. We required that a shower does not contain 
any "hot" crystals, is not matched to the projection of a charged track into the CC, 
and has E9/E25 OK. E9/E25 is the ratio between the total energy from 9 central 
crystals around the highest energy crystal and the total energy from the 25 crystals 
centering at the highest energy crystal in the calorimeter. E9/E25 OK means that 
for a shower with a certain energy, the ratio is within 99% of the E9/E25 photon hke 
distribution at this energy. 

In selecting candidate transition photons for the decay ip{2S) — > ■yrjc{2S), we 
required the energy to be in the range — [30 — 110] MeV. Because of the 
degraded photon resolution and efficiency in the endcap, we selected showers in the 
good barrel only, | cos^| < 0.81. 

3.3.3 Selection 

We required the Kg flight significance to be greater than 3. Kg flight significance is 
the ratio of the distance between the position of Kg where it decays to two charged 
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tracks and the beam spot to the measurement error of the position determined by 
the error matrix. In another word, the intersection of the two charged tracks that 
are Kg daughters was required to be far enough from the interaction point. 

The mass of a Kg candidates was required to be within 10 MeV of the Kg nominal 
mass, i.e., |M^o - 0.498 GeV| < 10 MeV. 

3.3.4 Neutral Pion Selection 

For selecting photons used to compose a 7r° or 77, CLEO-c software PhotonDecaysProd 
was used with the following defaults for the CLEO-c software packages NavPiOToGG 
and NavEtaToGG: E9/E25 OK, unmatched shower in good barrel or good endcap, 
and > 30 MeV. 

tt'^ candidates were obtained from the table of tt^ candidates that was produced 
by NavPiOtoGG. The default photon energy threshold of NavPiOToGG table is 30 
MeV in the good barrel and 50 MeV elsewhere. No hot crystals and E9/E25 OK 
were also required. The only additional requirement for selecting 7r°'s was to require 
the pull mass of the to be within ±3a. The pull mass is defined as the ratio of 
the difference between the mass of the 7r° candidate and the 7r° nominal mass divided 
by the standard deviation of the 7r° mass distribution. This means that the mass of 
the 7r° candidate was required to be close enough to the tt" nominal mass. 

3.3.5 r] Selection 

The 77 candidates were reconstructed in two ways, 77 — > 77 and 77 — > 7r"'"7r~7r°. 
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The ?7 — i> 77 candidates were obtained from the NavEtaToGG (CLEO-c software 
package) table similar to selecting the tt^ candidates. The default requirements of 
photons in NavEtaToGG are the same as NavPiOToGG. Similarly, the only additional 
requirement for selecting — 77 decays was to require the pull mass of the rj to be 
within ±3cr. 

For T] TT+TT^Tr", we required the mass of tt+tt^tt^ combination to be within 
10 MeV of the nominal 77 mass, i.e., |M^+jr-7rO ~ ^r?| < 10 MeV. Figure [STT] shows 
simulated t] mass distributions before the cut was applied in the two rjci'^S) modes 

with Tj — i> TTTTTT^. 

3.3.6 T]' selection 

Among the seven modes we considered, only one mode has an rj' decay, ric{2S) —>■ vrvr?/', 
7]' 7r7rr7(77). For this mode we required the invariant mass of the hadrons that 
compose the r]' to be within 10 MeV of the expected rj' mass, \M{inTr]) — M{r]')\ < 
10 MeV. Figure 13.21 shows simulated 77' mass distributions for this mode before the 
cut was applied. 

3.3.7 Suppression of Background Charmonium States 

In following section, we describe selection criteria used to suppress production of other 
charmonium states observed in the generic ijj{2S) MC samples. 
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Figure 3.1: MC simulations of r] candidate masses in ip{2S) events with the final 
states with Tnrr]{innr^) (top) and KKr]{7nn:^) (bottom). The solid histogram is the 
5 times luminosity background MC samples, while the dashed histogram is signal 
MC, arbitrarily scaled for clarity. The arrows show the selection cuts that were 
applied. All other event selection criteria have been applied. 
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Figure 3.2: MC simulations of r]' candidate masses in ^{2S) events with the final 
state vrvrr/', r/' — > 7r-7r7/(77). The solid histogram is the 5 times luminosity background 
MC samples, while the dashed histogram is signal MC, arbitrarily scaled for clarity. 
The arrows show the selection cuts that were applied. All other event selection 
criteria have been applied. 
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Hadronic Invariant Mass Difference 

To suppress decays of direct ^^{23) and Xc2, produced via ip{2S) — * 7Xc2, to tlie 
liadronic final states under investigation, we require tfie liadronic invariant mass dif- 
ference AM = M{iIj{2S))-M{X) to be in tlie range [0,100] MeV, wliere M(V^(2S)) = 
3686.1 MeV is tlie PDG value of the ip{2S) mass [7], and M(X) is the reconstructed 
hadronic invariant mass. 

7] Recoil Mass 

In selecting rjc{2S) decays to final states with an i] decay, large backgrounds arise 
from ip{2S) J/il){lS)ri decays. The branching fraction of this decay mode is 
(3.09 ± 0.08)% [7J. Therefore we applied t] recoil mass cuts to the modes rjc{2S) — >• 
7r7r?7(77) and 7nTri[innT^). In the mode vrvrr/', r]' 717177(77), although an rj is in the 
decay, we did not apply the t] recoil mass cut because the selection of r]' has already 
eliminated i'{2S) J/ipi^lS)!]. 

The 7] recoil mass cut was mode dependent. For mode rjc{2S) 7171^7(77), we 
rejected events with the rj recoil mass within 40 MeV of the J/ip mass by requiring 
\RecM{r]) - M{J/ip)\ > 40 MeV. For mode 'q^{2S) 7r7r77(7r7r7r°), we required 
|RecM(77) — M{J/'il))\ > 20 MeV. Simulated t] recoil mass distributions before 
applying this cut are shown in Figure 13. 3[ These cuts were found to be unnecessary 
for the T]c{2S) — > KKrji^'-)) and KKrjiT^Tnx^) modes. 
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Figure 3.3: MC simulations of the rj recoil mass distributions for ip{2S) events with 
final states -^"7^7(77) (top) and 7r7rr/(7r7r7r°) (bottom). The solid histogram is the 
5 times luminosity background MC samples, while the dashed histogram is signal 
MC, arbitrarily scaled for clarity. The arrows show the selection cuts that were 
applied. All other event selection criteria have been applied. 
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Photon Recoil Mass 

In modes that contain rj ^ 'j'y and vr" 77 decays, a large amount of background 
comes from using the photon from the decay ip{'2S) —>■ 'JXcJ with another photon. 
By applying mode dependent recoil mass cuts, we effectively remove a large amount 
of this background. 

For modes ric{2S) -mirj{^'~^) and rjc{2S) KKrj^'yy), since M{ri) = (547.51 ± 
0.18) MeV [7], a relatively large mass compared to the energy of the transition photons 
in the ipi^S) IXcJ decays (100-300 MeV), the transition photon is sometimes used 
as the low energy photon in the rj decay. Therefore we examined the photon recoil 
mass for the lower energy photon in 77 — 77 decays. For these two modes, we rejected 
events by requiring 

• |RecM(?] low energy 7) - M{xc2)\ > 20 MeV 

• \RecM{r] low energy 7) - M{xci)\ > 20 MeV 

• |RecM(r/ low energy 7) - M{xco)\ > 30 MeV 

For the mode 7/^(25) ^ KKti'^, since M(7rO) = (134.9766 ± 0.0006) MeV [7], a 
relatively small mass compared to the energy of the transition photon (100-300 MeV), 
the high energy photon is more likely swapped with the transition photon. Therefore 
we examined the photon recoil mass for the higher energy photon in 7r° 77 decays. 
For this mode, we rejected events by requiring 

• |RecM(7r° high energy 7) - M{xc2)\ > 20 MeV 

• |RecM(7r° high energy 7) - M(xci)| > 20 MeV 
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• |RecM(7r° high energy 7) - M{xco)\ > 30 MeV 

Similar to the case of the t] recoil mass cut, modes ric{2S) TnrrilTnnT^) and 
KKrilnmr^), even though they contain 7r° decays, the photons used to form the 7r° 
candidate did not contain contamination from ip{2S) 'JXcJ decays. Therefore no 
photon recoil mass cut was applied to these two modes. 

Figures [33] and [33] show the simulated distributions of photon recoil mass for the 
decay modes imri^'j'y), KKrj^'yj), and KKn^ before applying these cuts. 

J/-ip Rejection - tttt Recoil Mass 

For the modes that include charged pion pairs, we rejected events which contained 
decays of (25) n^n^J/ip. Specifically, we required |RecMjr+7r-— > 20 MeV 
for modes Gvr, KKtctt, KKmrir^, KKAn, and KsKStt, and RecM^+^- — Mj/^ < 
—30 MeV for mode An. With these requirements, we can suppress background from 
the decay ipi^S) n^n^J/tp, which is a dominant decay mode of 4'{2S) {B ~ 31.8%). 

The distributions of vrvr recoil mass for modes 47r, Gtt, KKtitt, KKmnT^, KKAtc, 
KsKSn before applying this cut are shown in Figures 13.61 through 13. 8[ 

J/ip Rejection - Hadron Invariant Mass 

Even after we rejected events with the tttt recoil mass cuts, J/ip —>■ {X — 2tt) decays, 
where X is a particular final state under study, were still observed in the generic 
ip{2S) MC sample. This background can be further suppressed by requiring the 
invariant mass of hadrons (excluding the tttt pair) to be displaced from the J/ if) mass. 
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Figure 3.4: MC simulations of distributions of low energy photon recoil mass for 
il){2S) events with final states -Kirrji^^) (top) and KKr]{'y'y) (bottom). The solid 
histogram is the 5 times luminosity background MC samples, while the dashed 
histogram is signal MC, arbitrarily scaled for clarity. The arrows show the selection 
cuts that were applied. All other event selection criteria have been applied. 
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Figure 3.5: MC simulations of distributions of low energy photon recoil mass for 
ip{2S) events with the final state KKtt^. The solid histogram is the 5 times lumi- 
nosity background MC samples, while the dashed histogram is signal MC, arbitrarily 
scaled for clarity. The arrows show the selection cuts that were applied. All other 
event selection criteria have been applied. 
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Figure 3.6: MC simulations of distributions of vrvr recoil mass for ip{2S) events with 
final states 47r (top) and Gvr (bottom). The solid histogram is the 5 times luminosity 
background MC samples, while the dashed histogram is signal MC, arbitrarily scaled 
for clarity. The arrows show the selection cuts that were applied. All other event 
selection criteria have been applied. 
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Figure 3.7: MC simulations of distributions of tttt recoil mass for ip{2S) events 
with final states KKtttt (top) and KKTrinr^ (bottom). The solid histogram is the 
5 times luminosity background MC samples, while the dashed histogram is signal 
MC, arbitrarily scaled for clarity. The arrows show the selection cuts that were 
applied. All other event selection criteria have been applied. 
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Figure 3.8: MC simulations of distributions of tttt recoil mass for tp{2S) events 
with final states KK4tt (top) and KsKStt (bottom). The solid histogram is the 
5 times luminosity background MC samples, while the dashed histogram is signal 
MC, arbitrarily scaled for clarity. The arrows show the selection cuts that were 
applied. All other event selection criteria have been applied. 
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More specifically, for the modes X = Gvr, KKtttt, KKtttt^, KKAtt, KsKStt, we reject 
events by requiring \M{X - 2n) - M{J/ij)\ > 30 MeV. 

Figure [33] shows the distributions of hadron invariant mass for the 6n mode before 
applying the cut. The distributions for other four modes {KKtttt, KKt:tttt^, KKAtt, 
and KsKStt) before applying the cut are shown in Figure I3.10[ 

3.3.8 Photon Background Suppression 

We investigated background suppression strategies related to the relationship between 
the candidate transition photon and the charged tracks in the event. 

Photons which are emitted from charged particles are typically in the same direc- 
tion as the initial momentum vector of the charged particle at the IP. These photons 
are called final state radiation (FSR). If a shower is close to a nearby track in the 
CC, it may be a splitoff shower (a shower that split off from a hadronic interaction 
of a nearby track in the CC). To identify and remove these backgrounds, the angle 
between the transition photon and the initial momentum of the closest pion (cos6'^,7r) 
and the distance between the transition photon and the nearest track (d^^trk) were 
studied. Although either of these can help remove FSR and splitoff showers, a cut 
on the angle between the transition photon and the initial momentum of the closest 
pion is more efficient in suppressing final state radiation and a cut on the distance 
between the transition photon and the nearest track is more efficient in removing 
splitoff showers. In this analysis, the following rules were applied while making the 
decision on which of these two cuts would be applied to a mode: 
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Figure 3.9: MC simulations of distributions of hadron invariant mass M{X — 27r) for 
ip{2S) events with the final state 6tt. The solid histogram is the 5 times luminosity 
background MC samples, while the dashed histogram is signal MC, arbitrarily scaled 
for clarity. The arrows show the selection cuts that were applied. All other event 
selection criteria have been applied. 
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Figure 3.10: MC simulations of distributions of hadron invariant mass M{X — 2ti) 
for ip{2S) events with final states KKtttt (top left), KKninT^ (top right), KKAtt 
(bottom left) and KsKStt (bottom right). The solid histogram is the 5 times lumi- 
nosity background MC samples, while the dashed histogram is signal MC, arbitrarily 
scaled for clarity. The arrows show the selection cuts that were applied. All other 
event selection criteria have been applied. 
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• If the hadronic decay only contained charged pions, cut on the angle between 
the initial pion momentum and the candidate transition photon. 

• If the hadronic decay only contained charged kaons, cut on the distance of the 
closest track to the candidate transition photon. 

• If the hadronic decay contained both charged pions and kaons, cut on the dis- 
tance of the closest track to the candidate transition photon, then we examined 
if a cut on the angle between initial pion momentum and the candidate transi- 
tion photon improves the S"^ / {S + B). (We did not apply both types of cuts to 
the same mode.) 

We optimized the selection criteria by evaluating the figure of merit S^/ {S + B) 
at different cut values. The value of S was from the signal MC sample with the 
assumption that B{^{2S) -fVc{2S)) x B{ri^{2S) X) = 2.6 x lO'^ for all ri^{2S) 
decay modes. The value of B was from the sum of the 5 times generic ip{2S) and 
continuum MC samples scaled to the true luminosity of data. The plot of S"^/ (S + B) 
versus cut value was examined to determine the optimal cut value. An example of 
the S^/{S + B) is shown in Figure 13.111 for the Att mode. The S'^/{S + B) for the 
other modes are shown in Figures 13.121 through I3.14[ The optimized cut for each 
mode is listed in Table 13. 2[ 
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Figure 3.11: MC study of S'^/{S + B) for different cuts on tlie angle between the 
transition photon candidate and the closest pion track for the An mode. The arrow 
shows the cut value that was selected. 
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Figure 3.12: MC study of 5^/(5 + B) for different cuts on the distance between a 
candidate transition photon shower and the nearest track in the CC or the angle 
between the transition photon candidate and the closest pion track for modes Qtt 
(top left), KKiTTT (top right), KKtt^ (bottom left), and KsKtt (bottom right). The 
arrows show the cut values that were selected. 
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Figure 3.13: MC study of 5^/(5 + B) for different cuts on the distance between a 
candidate transition photon shower and the nearest track projection into the CC 
or the angle between the transition photon candidate and the closest pion track 
for modes Tnrr]{jj) (top left), mn]{7n:Tr^) (top right), KKr]{'y^) (bottom left), and 
KKr](TnnT^) (bottom right). The arrow shows the cut value that was selected. 
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Figure 3.14: MC study of S'^/(S + B) for different cuts on the distance between a 
candidate transition photon shower and the nearest track projection into the CC 
or the angle between the transition photon candidate and the closest pion track for 
modes vrvrr/' with rj' —>■ TTirrj^-jj) (top left), KKinnr^ (top right), KKAtt (bottom 
left), and KsKStt (bottom right). 
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3.3.9 Global Event Selection 
4-Constraint Kinematic Fit x^/dof 

To select events consistent with full detection of ip{2S) — > ■yric{2S),ric{2S) ^ X we 
compute the 4-constraint kinematic fit x^/dof using the package FitEvt [80j . Both 
momentum and vertex constraints were imposed. For each mode, a vertex fit of the 
charged tracks from the IP was performed first and the fit result was used as the 
origin of the photon(s) in the kinematic fit. In the fit we imposed the constraint 
that the total 4-momentum of all tracks and neutral particles be equal to the total 
4-momentum of the ip{2S). We optimized a cut on the resulting overall x^/dof in 
the same manner as the optimization of the cuts on nearby tracks (see Section [3. 3. 8p . 
using S"^/ {S + B) as the figure of merit. An example of the 5*^/ (5* + B) optimization 
for the 4-momentum kinematic fit x^/dof is shown in Figure [3.151 for the 47r mode. 
The S'^/{S + B) for the other modes are shown in Figures [3.161 through 13.181 The 
optimized cuts are listed in Table 13.21 

4-Constraint Kinematic Fit x^/dof of Hadrons Only 

Similar to the 4-constraint kinematic fit for the whole event, we used the FitEvt 
package [80] to perform the fit only for the hadrons of 4'{2S) X candidates, 
both momentum and vertex. In the fit we impose the constraint that the total 
4-momentum of all hadrons be equal to the total 4-momentum of the ip{2S). If a 
candidate is a ip{2S) — > '-fric{2S) event, then the 4-C kinematic fit x^/dof of the whole 
event should peak at zero, and that of just the hadrons should be distributed away 
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Table 3.2: Optimization of 4-C kinematic fit x^/dof, along with the transition pho- 
ton distance to the nearest track (d^^trk) or the angle between the transition photon 
and the closest pion (cos^-y^Tr). 



Mode 


d-y^trk (cm) or cos Ory^^^ 




47r 


COs6'-y,:;r < 0.93 


< 4.5 


Gtt 


- 


< 5.0 


KKtttt 


d'y,trk > 45 


< 4.0 




dj,trk ^ 35 


< 4.0 


KsKtt 




< 3.5 


777777(77) 




< 2.0 


777777(777777°) 


COS ^-y,^ < 0.93 


< 3.0 


777777', v' ~^ 777777(77) 




< 3.0 






< 3.5 






< 5.0 






< 2.5 


KKAn 




< 4.0 


KsK3tt 




< 4.0 



98 



471 




123456789 10 
4C fit x^/dof 



Figure 3.15: MC study of S'^ /(S+B) for different cuts on the 4-momentum kinematic 
fit x^/dof for the Air mode. The arrow shows the cut value that was selected. 
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Figure 3.16: MC study of S'^/ (S+B) for different cuts on the 4-momentum kinematic 
fit x^/dof for modes Gvr (top left), KKmr (top right), KKn^ (bottom left), and 
KsKiT (bottom right). The arrows show the cut values that were selected. 
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Figure 3.17: MC study of 5^/ (S+B) for different cuts on the 4-momentum kinematic 
fit x^/dof for modes 7r7rr/(77) (top left), 7nTr]{TnnT^) (top right), KKr]{'-f'-f) (bottom 
left), and KKr]{7nnr^) (bottom right). The arrows show the cut values that were 
selected. 
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Figure 3.18: MC study of 5^/ (5+5) for different cuts on the 4-momentum kinematic 
fit x^/dof for modes vrvrr/' with ry' — > 7r7rr/(77) (top left), KKti'kti^ (top right), i^i^47r 
(bottom left), and KsKSn (bottom right). The arrows show the cut values that 
were selected. 
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from zero. This cut was evaluated, but was found not to suppress the background as 
efficiently as hoped. Therefore it was excluded from the final set of selection criteria. 

3.4 Signal Extraction Procedures 

In this section, we describe the procedure we developed for extracting the signal 
yields and resonance parameters from the 25.9 M ip{2S) decays. We fit the mea- 
sured, or unconstrained, photon energy distribution in the range E^^ = [30, 94] MeV. 
Figure 13.191 shows the expected signal distribution for the modes Att and KsKtt 
for a mock data sample of 25.9 M ip{2S) decays generated under the assumption 
that the partial widths for ric{2S) X are the same as for ric{lS) — > X, and 
B{iIj{2S) 7r7e) = 2.6 x 10"^ 

3.4.1 Detector Resolution 

Since we extract the signal from a measured photon distribution, a Crystal Ball 
function [HTl [82] is used to fit the detector resolution distribution from the signal 
MC samples described in Section I3.2.3[ The Crystal Ball function (named for the 
Crystal Ball experiment) is designed to provide a good representation of the photon 
energy distribution as measured with a high-precision crystal calorimeter like CLEO- 
c's. The resolution is defined as the standard deviation of the difference between the 
measured and generator-level photon energy. The RooFit software package [83j was 
used to perform these fits. 
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Figure 3.19: The expected signal distributions for modes An (top) and KsKtt (bot- 
tom) for a simulated data-sized sample based on the assumed product branching 
fractions listed in Table [L5l The solid histogram is the sum of the 10 times generic 
ip{2S) and 5 times continuum MC samples scaled to one time luminosity of data. The 
points are randomly thrown signal MC events, normalized to the product branching 
fraction expectation listed in Table II. 5i The signal MC was generated with the 
expected signal mean of 48 MeV and width of 14 MeV. 
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The Crystal Ball function is defined as 




) 



—n 



for 



for < —a 



(7 



> —a 



(3.1) 



where is the normalization factor and a, n, x, and a are parameters that govern 
the shape. The parameter x represents the mean of the resolution function and a is 
the resolution in the photon energy. The Crystal Ball function features a Gaussian 
distribution with a power-law tail (index n) at the low end below a threshold (given 
by the parameter a) to account for energy leakage within the calorimeter. 

To obtain successful fits for the resolution function, the parameter n was fixed at 
140 for all modes. It was observed that when n varied over the range from 40 to 140, 
the shape of the Crystal Ball function and the fit showed very little variation. 

As an example Figure 13.201 shows the resolution fit of the Air mode. Resolution 
fits of other modes are shown in Appendix \M 

3.4.2 Background Parametrization 

A necessary ingredient for obtaining reliable results was the use of an appropriate 
function to describe the background shape. For this purpose we can either use a 
polynomial or other function to parametrize the background or rely on a MC sim- 
ulation that at least approximates the underlying physics. Clearly the latter is 
potentially more reliable, but it is first necessary to study the reliability of our MC 
samples and the composition of the background. An independent data sample is 
required for such a study, and we used ^'(25') n^n^J/ip events with four specific 
J/ip decay modes. These modes are J/ip ^ An, KKtttt, KKtt^, and KsKtt. Note 
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Figure 3.20: Resolution fit of decay il}{2S) 7??c(25'), r]c{2S) An. 
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that these modes gave us a sample of decays modes with vr, K, 7r°, and Ks- When 
selecting J/ip events, the vrvr recoil mass was required to be within 20 MeV of the 
mass of J/ip. All other event selection criteria for ip{2S) '~fr]c{2S) decays were 
applied, except that the vrvr recoil mass and J/ip ^ {X — 2ti) rejection cuts were 
removed and the hadronic invariant mass cut was adjusted for the J/ip mass instead 
of the mass of 'il){2S). 

The two primary sources of low energy shower backgrounds are splitoff showers and 
FSR. The splitoff showers are generally very close to the track associated with them. 
Most splitoff showers are removed by the cut on the transition photon candidate 
distance to the nearest track, but some pass through the cut. Additional splitoff 
showers can be removed with a tighter cut on the full event 4-constraint kinematic fit 
X^/ dof . Figure 13^211 shows the fits of the measured photon energy distribution for data 
in J/ip — i> 47r decays with an accompanying (background) photon candidate. The 
upper left plot in the figure shows a fit using the 5 times luminosity generic ip{2S) and 
continuum MC samples with one free parameter (normalization) for the background. 
The upper right plot shows a two parameter fit to the separate distributions for 
non-splitoff showers (mostly FSR) and for splitoff showers. A shower was MC-truth 
tagged as a splitoff shower if its parent was not a photon. ("MC-truth tagging" 
uses detailed information stored at the time the MC was created to select only those 
particles that come from a particular process.) The lower left plot shows a histogram 
(one parameter) fit for events when the 4-constraint kinematic fit x^/dof has been 
tightened to be less than 3. The lower right plot shows a fit with a linear background 
function. Figures [3.221 through 13.241 show similar plots for the other three J/ip decay 
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Figure 3.21: Background fits of the measured photon-energy distribution for 
ip{2S) — > vr^vr^J/'i/^ events with J/ij: — > Att and a low energy shower. The up- 
per left plot is the fit of data to a single MC histogram. The upper right plot is 
the fit of data to two MC histograms with independent parameters, one (dashed) 
histogram excluding splitoff showers and the other (dotted) including only splitoff 
showers. The lower left plot is the fit of data to a single MC histogram with a 
tighter cut of full event 4-C fit x^/dof < 3.0. The lower right plot is the fit of data 
with a linear function. 
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Figure 3.22: Background fits of tlie measured photon-energy distribution for 
ip{2S) — > 7r^7r~J/'0 events with J/ip ^ KKtttt and a low energy shower. The 
upper left plot is the fit of data to a single MC histogram. The upper right plot 
is the fit of data to two MC histograms with independent parameters, one (dashed) 
histogram excluding splitoff showers and the other (dotted) including only splitoff 
showers. The lower left plot is the fit of data to a single MC histogram with a 
tighter cut of full event 4-C fit x^/dof < 3.0. The lower right plot is the fit of data 
with a linear function. 
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Figure 3.23: Background fits of tlie measured photon-energy distribution for 
ip{2S) — > TT^TT^J/ip events witli J/^p — > KKtt^ and a low energy shower. The 
upper left plot is the fit of data to a single MC histogram. The upper right plot 
is the fit of data to two MC histograms with independent parameters, one (dashed) 
histogram excluding splitoff showers and the other (dotted) including only splitoff 
showers. The lower left plot is the fit of data to a single MC histogram with a 
tighter cut of full event 4-C fit x'^/dof < 3.0. The lower right plot is the fit of data 
with a linear function. 
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Figure 3.24: Background fits of tlie measured photon-energy distribution for 
^{2S) — > TT^TT^J/tp events with J/ip ^ KsKtt and a low energy shower. The 
upper left plot is the fit of data to a single MC histogram. The upper right plot 
is the fit of data to two MC histograms with independent parameters, one (dashed) 
histogram excluding splitoff showers and the other (dotted) including only splitoff 
showers. The lower left plot is the fit of data to a single MC histogram with a 
tighter cut of full event 4-C fit x^/dof < 3.0. The lower right plot is the fit of data 
with a linear function. 
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modes. The fits show that the low energy shower distribution of the background MC 
sample represents the data quite well. Fitting with the two sources of backgrounds 
separately does not lead to obvious improvement of background fitting based on the 
resulting x^- Table [373] shows the results of the four types of fits for all four J/ip 
decay modes studied. 

Table 3.3: Results of x^/dof for background fits of the J /ijj decay modes. In this 
table, "1 hist" means the results of fitting data with MC histogram; "2 hist" means 
the results of fitting data with two MC histograms, one for non-splitoff histogram 
and one for splitofF histogram; "tight x^/dof" means the results of fitting data with 
MC histogram when a tighter 4-C fit x^/dof cut was applied; "linear" means the 
results of fitting data with a linear function. 



Mode 


xVdof 


1 hist 


2 hist 


tight 


linear 


An 


1.25 


1.31 


1.92 


1.27 


KKlTTT 


0.59 


0.61 


0.74 


0.99 


KKn^ 


0.59 


0.52 


0.96 


0.99 


KsKn 


1.99 


2.05 


1.79 


2.48 



This study showed that a linear background function provides a reasonable alter- 
native to the MC-predicted distributions but one that is not based on the physics 
of ipi'^S) decays or the CLEO-c detector. Based on this investigation, we chose 
the single-histogram as our standard background shape fit in studying the ^/'(2S') 
jrjc{2S), Tjci^S) — * X decays. The linear fit is included as an alternative for evalu- 
ating the systematic uncertainty in the fits. 
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3.4.3 Signal Fit 

The measured, or unconstrained, shower energy distribution was fit for a signal using 
a non-relativistic Breit-Wigner function (Cauchy distribution) [HHIHS] convolved with 
the Crystal Ball resolution function. The parameters of the Crystal Ball function, n, 
a, and a were obtained from the resolution fit and were fixed in the signal fits. The 
Breit-Wigner function is defined as follows: 

where F is the width of the resonance and x is the mean of the variable x. In our 
signal fits, x represents the measured photon energy. 

For the signal fits, the parameters are extracted with a binned hkelihood fit, 
minimizing the negative logarithm of hkelihood using RooFit [83]. For each mode 
the histogram background shape was the sum of events from the 10 times generic 
ip{2S) and 5 times continuum MC sample which passed our selection criteria. We fit 
the data with the above signal function and background histogram. 

To verify the fitting procedure, the decays of i'{2S) — > 7Xc2, Xc2 X were studied. 
The details of this study can be found in Section 13.51 

3.5 Test of Fit Procedure with '0(25') lXc2 Decays 

The study of the decay channel 'ipi^S) 7Xc2 was used to verify the analysis proce- 
dures for the decay ip{2S) — > 'yric{2S). 
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3.5.1 Introduction 

The radiative decay of ipi^S) — ^ 7Xc2 is similar to the decay ipi^S) 'yric{2S). The 
Xc2 is the P-wave spin triplet charmonium state with total angular momentum J = 2 
mass 3556.20 ±0.09 MeV, and width 2.05 ±0.12 MeV ^Tj. For a iIj{2S) meson nearly 
at rest, the transition photon from the decay ip{2S) IXdi should have energy 
around 127.60 MeV. This is 80 MeV higher than the energy of the transition photon 
in the decay '^{2S^ — > '-)ric{2S). While the photon energy range and background 
parametrization for the Xc2 study are different from the rjc{2S) study, the other event 
features and analysis procedures are very similar. 

3.5.2 Analysis Procedures 

Since the purpose of studying this decay channel is to evaluate the analysis procedure 
of V'(2>S') — i> ■yric{2S), similar procedures were applied to 4'{2S) lXc2 decays. The 
results were obtained from similar resolution and signal fits on signal MC and data 
samples. 

Same as t]c{2S) decays, the data sample for this study was the full sample of Data 
32 and Data 42 of CLEO-c data, 25.9 M tlj{2S) decays. As described in SectionE231 
the Xc2 signal MC sample was generated to replicate the full CLEO-c ^^(25*) sample, 
similar to the ric{2S) sample. 

To validate the event selection criteria for il){2S) —>■ 'yric{2S) decays, the same 
criteria for 4'{2S) 7?7c(2S') decays were applied to 4'{2S) lXc2 decays, except 
that AM = M^(2s) - Mhad = (0, 100) MeV was removed and = (30, 110) MeV 
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was replaced with = (90, 145) MeV. 
3.5.3 Results 

For each mode, the Crystal Ball function [82j was used to determine the detector 
resolution of the measured photon energy from signal MC. The parameter n of the 
Crystal Ball function was set at a fixed value of 140 in order to obtain successful fits. 
From the resolution fits, the value of a and a were determined. These values were 
then fixed and applied in fitting the measured photon energy distribution of data. 

In the signal fit of each mode, we used a Breit-Wigner distribution convolved 
with a Crystal Ball resolution as the signal function. In addition to fixed n, a 
and cr, the width (F) was fixed to 2.05 MeV. Unlike the ip{2S) — > ■yric{2S) study, the 
background shape was either a linear function, constant, or no background, depending 
on the fitting status of the linear background fit. For all modes, a fit with a linear 
background and mean left free was performed first. If the linear fit did not converge, 
and the signal mean was within ±1 MeV of the nominal mean (127.60 MeV), then a 
fit with a constant background and mean left free was attempted. If the linear fit 
did not converge and the mean was not within ±1 MeV of the nominal value, then 
a fit with a constant background and the mean fixed to 127.60 MeV was attempted. 
If the linear fit resulted in a background with a negative yield, then a fit with no 
background was performed. According to these rules, the modes 4it, Gtt, KKtctt, 
KKAtt and KsKStt were fit with a linear background; the modes KKtt^, 7T7n]{'j'-f) , 
TnTrilTTTT-K^) and KKtctitt^ were fit with a constant background; the modes KKri{'y'y) 
and TTTTT]' with rj' nnri^'y'y) were fit with a constant background and fixed mean; 
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and the modes KsKir and KKriiji-KTi^) were fit without a background. 

The efficiencies and the resolution fit results are listed in Table [3^ The branching 
fractions of the decays Xc2 X are given by 

■ 

B{Xc2 - X) = — ^ (3.3) 

where N^ig is the number of signal events (yields) determined by the signal fits, e is 
the efficiency, N^(2S) is the total number of ipi^S) decays, and B{4'{2S) 1X02) is 
the branching fraction of ip{2S) 7Xc2- When calculating the branching fractions 
of Xc2 X, we used A^^(2S) = 25.9 x 10^ p] and B{tp{2S) ^ 7x^2) = (8.1 ± 0.4)% 
[7j. The results of the signal fits and the corresponding ip{2S) 7Xc2 branching 
fractions are listed in Table 13.51 The data fits for all thirteen modes are shown in 
Figures 13.251 through 13.311 while the resolution fits are shown in Appendix |X1 

3.5.4 Discussion and Conclusion for the Xc2 Study 

From the results we can see that for most modes the results are consistent with the 
world averages. However, two modes show inconsistency with the PDG values, Gtt 
and 71717]' . 

The previous CLEO measurement [871 EH] shows the result of the branching frac- 
tions of Xc2 X decays to be (0.49 ± 0.12 ± 0.05 ± 0.03) x lO^^ for Trvrr/ mode, 
< 0.33 X 10^3 for KKt] mode, (0.51 ± 0.18 ± 0.05 ± 0.03) x 10"^ for Trvrr]' mode, and 
(0.31 ± 0.07 ± 0.03 ± 0.02) x 10"^ for KKt^^ mode. Our results are consistent with 
this measurement except for Xc2 ttttt]' mode. The cuts in our measurements are 
not optimized for Xc2 studies. 
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Table 3.4: Efficiencies and resolution function parameters for ■0(25') — > 7Xc2 decays. 



Mode 


e{%) 


Res a (McV) 


XBall a 


An 


26.49 ±0.17 


6.23 ±0.04 


1.26 ±0.02 


Gtt 


18.43 ±0.15 


6.40 ±0.05 


1.33 ±0.03 


KKttt: 


24.91 ±0.16 


6.30 ± 0.04 


1.42 ± 0.03 


KK7r° 


25.48 ±0.16 


6.12 ±0.05 


1.10 ±0.02 


KsKtt 


25.46 ±0.16 


6.38 ± 0.04 


1.56 ±0.04 


777777(77) 


16.68 ±0.14 


5.97 ±0.05 


1.33 ±0.04 




11.77±0.12 


6.24 ±0.06 


1 40+0-05 

-'-•^"-'-0.04 


TTTrrj' 


12.86 ±0.15 


6.13 ±0.08 


1 26+°°^ 

-■-■^"-0.04 


KKrj{jj) 


18.86 ±0.15 


6.12 ±0.05 


1.16 ±0.03 




14.01 ±0.13 


6.06 ±0.06 


1.05 ±0.03 




11.93 ±0.12 


6.28 ±0.06 


1 fiq+0.08 


KKAtt 


12.97 ±0.13 


6.37 ±0.06 


1.45 ± 0.05 




11.21 ±0.13 


6.37 ± 0.06 


1 in+0.0.5 
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Table 3.5: Results for '0(25') jXc2 study. The value B{iIj{2S) 7X02) = 
(8.1 lb 0.4)% was used to determine B = B{xc2 X). Only statistical errors are 
listed. Note that the PDG value for KsKtt listed here is corrected, because the 
PDG did not convert KsKn to K^Kn when listing BES-II result [3 [86]. The mean 
in the 7r7rTy(7r7r7rO) and Trvrry' fits was fixed to 127.60 MeV, and all others are free. 



Mode 


Mean (MeV) 


'Slg 


B (10-3) 


BpBG (10-3) 


4:71 


127.23 ±0.10 


7215 ± 119 


13.0 ±0.2 


12.5 ± 1.6 


Qn 


127.09 ±0.11 




15.7 ±0.3 


8.7± 1.8 


KKtctc 


127.20 ±0.12 


4717 ±95 


9.0 ±0.2 


10.0 ±2.6 




127.71 ±0.58 


219 ± 17 


0.41 ±0.03 


0.36 ±0.09 


KsKtt 


127.87 ±0.43 


294 ± 17 


0.80 ±0.05 


0.71 ±0.11 


7r7rr](77) 


128.15l°;f, 


97± 12 


0.70 ±0.10 


0.56 ±0.15 




126.86+2-10 


31 ±7 


0.55 ±0.13 


0.56 ±0.15 


rnirj' 


127.60 


3.7 ±5.2 


0.08 ±0.11 


0.59 ±0.22 




127.60 


29 ±8 


0.19 ±0.05 


< 0.4 


KKri{'K'K'K^) 


127.96l|2} 


17±5 


0.26 ±0.08 


< 0.4 


KKmrTT^ 


127.39 ±0.14 


3197 ±62 


12.8 ±0.3 




KKAti 


127.24 ±0.18 


2249 ± 68 


8.3 ±0.3 






127.49 ±0.23 


14531^^ 


7.2 ±0.3 
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Signal Fit - 47i 




Figure 3.25: Measured photon energy for the decay mode ^^{23) — > jXc2, Xc2 — > 47r. 
The points are from the 25.9 M ip{2S) data sample. The dashed hne is the result 
of the signal fit, the dotted line is the background fit, and the solid line is the sum 
of the signal and background fits. 
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Signal Fit - 671 




Signal Fit - KKtiti 




Figure 3.26: Measured photon energy for the decay modes tp{2S) 7Xc2iXc2 — *■ Gvr 
(top) and KKtttt (bottom). The points are from the 25.9 M ip{2S) data sample. 
The dashed hne is the result of the signal fit, the dotted line is the background fit, 
and the solid line is the sum of the signal and background fits. 
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Figure 3.27: Measured photon energy for the decay modes ■0(25) 7Xc2iXc2 
KKtt^ (top) and KsKtt (bottom). The points are from the 25.9 M tp{2S) data 
sample. The dashed Une is the result of the signal fit, the dotted line is the 
background fit, and the solid line is the sum of the signal and background fits. 
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Figure 3.28: Measured photon energy for the decay modes ip{2S) 'yXc2,Xc2 
7r7rry(77) (top) and Tnrr](Tnnr^) (bottom). The points are from the 25.9 M tp{2S) 
data sample. The dashed hne is the result of the signal fit, the dotted line is the 
background fit, and the solid line is the sum of the signal and background fits. 
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Figure 3.29: Measured photon energy for the decay modes ip{2S) lXc2-,Xc2 
KKr]{'~f'y) (top) and KKT]{'innT^) (bottom). The points are from the 25.9 M tpi'^S) 
data sample. The dashed hne is the result of the signal fit, the dotted line is the 
background fit, and the solid line is the sum of the signal and background fits. 
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Figure 3.30: Measured photon energy for the decay modes ip{2S) — > 7Xc2)Xc2 ~^ 
vrvrr/', -q' TTm]{'^^) (top) and KKninT^ (bottom). The points are from the 25.9 M 
tp{2S) data sample. The dashed hne is the result of the signal fit, the dotted line is 
the background fit, and the solid line is the sum of the signal and background fits. 
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Signal Fit - KK47i 
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Figure 3.31: Measured photon energy for the decay modes ip{2S) 7Xc2iXc2 
KKAtt (top) and KsK'i'K (bottom). The points are from the 25.9 M ?/;(25) data 
sample. The dashed hne is the result of the signal fit, the dotted line is the 
background fit, and the solid line is the sum of the signal and background fits. 
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The PDG branching fraction average for Xc2 Svr is dominated by measurement 
from the BES-I collaboration [89]. As can be seen in Table [3^ the BES-I results 
are systematically lower than our results, while our measurement of Xc2 —>■ 47r is 
completely consistent with the PDG average. 

Table 3.6: Detailed comparison of decay ip{2S) — > 7Xc2 results between our results, 
BES-I [89], and the PDG. The PDG branching fraction for Xc2 Air mode is from 
a 28 parameter fit using properties of the XcO, Xci, Xc2i V'(25') |7j. The BES-I 
measurement of Xc2 Stt dominates the PDG average. All branching fractions 
listed in units of 10^^. 



Mode 


PDG 


This result 


BES-I 


This Resuh/BES-I 


An 


12.5 ± 1.6 


13.00 ±0.22 


9.2 ±2.4 


1.4 ±0.4 


6n 


8.7± 1.8 


15.75 ±0.29 


8.7± 1.9 


1.8 ±0.4 


KKnn 


10.0 ±2.6 


9.03 ±0.18 


7.6 ± 1.9 


1.2 ±0.3 



For the mode nnrj' ^we were not able to explain for the discrepancy other than a 
statistical fluctuation in either our or the previous CLEO measurement [571 155]. 

Some of the decay modes like KKnnn^, KKAn and KsK3n have not been mea- 
sured before. For the decay mode KKt], only upper limits were given in PDG based 
on previous measurements. In this study, the branching fractions of Xc2 to these 
modes were determined. Even though the decays were studied with the optimized 
criteria for ip{2S) 'yric{2S) decays, the branching fractions results are consistent 
with values listed by the PDG. If the decays were studied with the optimized criteria 
for Xc2 decays, the efficiencies would have been higher and the statistical error would 
have been reduced. 
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Overall, the results of the study of 4'{2S) 7Xc2, Xc2 — ^ X decays establish the 
validity of the analysis procedures for ip{2S) 'yric{2S). There is no explanation for 
the small, but statistically significant, discrepancy observed in the vrvrr/' mode, but 
since the analysis of this mode is not appreciably different from the higher-statistics 
modes, we assume that it too has no fundamental flaw. 

3.6 Results for ^(25) ^ 7r]c{2S) 

The measured photon energy results for the rid^S) signal region are listed in Table [3?71 
The resolution function fits are shown in Appendix |X] and the corresponding values 
for a and a are fixed in the fits of the data sample. The ric{2S) resonance parameters 
M = 3638 MeV {E^ = 48 MeV) and T = 14 MeV are also fixed in the fits. The 
photon energy distributions for all eleven modes are shown in Figures 13.321 through 

EM 

Two methods are employed for determining the 90% confidence level upper limits 
on the yields listed in Table 13.71 For modes which do not contain t] decays, we 
have no apparent signals and substantial backgrounds. Therefore, we determine an 
upper limit by integrating the distribution defined from the nominal yield result (a 
bifurcated Gaussian) up to 90% of the area in the positive physical range. We use toy 
MC studies to verify that we get consistent results with the nominal fit integration. 
For modes with r] decays, a very limited number of events pass our selection criteria, 
either signal or background. Therefore, we use the Feldman and Cousins method 
|90] . The procedure of determining the upper limits is described in more detail in 
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Figure 3.32: Measured photon energy for the final states ^in (top) and jGtt (bottom) 
in the r]c{2S) signal region. The points are from the 25.9 M ip{2S) data sample. 
The dotted lines are a Breit-Wigner convoluted with the Crystal Ball resolution 
signal shape. The dashed lines are the background histogram determined from the 
10 times luminosity generic ip{2S) and 5 times luminosity continuum MC samples. 
The solid lines are the sum of the signal and background. 
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Figure 3.33: Measured photon energy for the final states '^KKtt^ (top) and jKsKtt 
(bottom) in the r]c{2S) signal region. The points are from the 25.9 M ip{2S) data 
sample. The dashed lines are a Breit-Wigner convoluted with the Crystal Ball 
resolution signal shape. The dotted lines are the background histogram determined 
from the 10 times luminosity generic ip{2S) and 5 times luminosity continuum MC 
samples. The solid lines are the sum of the signal and background. 
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Figure 3.34: Measured photon energy for the final states jKRittt (top) and 
'yKKTnnr^ (bottom) in the r]c{2S) signal region. The points are from the 25.9 M 
il>{2S) data sample. The dashed lines are a Breit-Wigner convoluted with the Crys- 
tal Ball resolution signal shape. The dotted lines are the background histogram 
determined from the 10 times luminosity generic ip{2S) and 5 times luminosity con- 
tinuum MC samples. The solid lines are the sum of the signal and background. 
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Figure 3.35: Measured photon energy for the final states ^KKAtt (top) and ^KsKStt 
(bottom) in the r]c{2S) signal region. The points are from the 25.9 M ip{2S) data 
sample. The dashed lines are a Breit-Wigner convoluted with the Crystal Ball 
resolution signal shape. The dotted lines are the background histogram determined 
from the 10 times luminosity generic '<p{2S) and 5 times luminosity continuum MC 
samples. The solid lines are the sum of the signal and background. 
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Figure 3.36: Measured photon energy for the final states jmrri (top) and ^KK-q 
(middle) and 77r7r77' (bottom) in the rjc{2S) signal region. The points are from the 
25.9 M ip{2S) data sample. The histograms are from the 10 times luminosity generic 
'0(25') and 5 times luminosity continuum MC samples. The solid red arrows enclose 
the signal region, while the dashed blue arrows enclose the sideband region. 
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Figure 3.37: Measured photon energy for the final state ^KKtt in the rjc{2S) signal 
region. The points are from the 25.9 M ip{2S) data sample. The dashed lines 
are a Breit-Wigner convoluted with the Crystal Ball resolution signal shape. The 
dotted lines are the background histogram determined from the 10 times luminosity 
generic ip{2S) and 5 times luminosity continuum MC samples. The solid line is the 
sum of the signal and background. 
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Appendix [Bl 



Table 3.7: Results for tp{2S) — > jrjc{2S), r]c{2S) X. The photon energy res- 
olutions were determined by fitting the resolution function from the rjc{2S) sig- 
nal MC sample. Other parameters were from signal fits of data. The statisti- 
cal significance, defined as V— 2Aln£, where is C is the likelihood, and the dif- 
ference is computed between the standard fit with both signal and background 
and an alternative background-only fit. The product Bi x B2 is defined as 
B{ip{2S) 7 r]c{2S)) X B{r]c{2S) X). All upper limits are at 90% confidence 
levels and only statistical errors are shown. 



Mode 


e 


Res (T 


XBall a 


Sgnf 






Bl X B2 


{X) 


(%) 


(MeV) 




i^) 


(Afob./Afbg) 




(10-6) 


A-K 


20.94 ±0.16 


4.68 ± 0.04 


1.50 ±0.04 


4.15 




< 64.8 


< 12.0 




14.71 ±0.14 


4.70 ±0.05 


1 43+0.05 
-^■^•-'-0.04 


0.56 


lo.iil?- 


< 36.6 


< 9.6 


KKn-K 


20.15 ±0.16 


4.57 ± 0.04 


1.45 ± 0.04 


0.82 


12.8tl^;8 


< 35.2 


< 6.7 


KK-K° 


18.82 ±0.15 


4.54 ± 0.04 


1 01 +0.04 
'^■■^'^-0.03 


1.44 


7 =+6.2 


< 16.0 


< 3.3 


KsK-K 


20.69 ±0.15 


4.61 ±0.04 


^■"^"-0.04 


0.94 


4.011° 


< 11.0 


< 3.0 


KKn 


7.89 ± 0.04* 


4.58 ± 0.03 


1.56 ±0.03 


1.73 




< 21.9 


< 10.7 


■KTTTj 


6.00 ±0.05* 








4 / 4.3 


< 4.3 


< 2.3 




8.63 ±0.11 








2 / 1.8 


< 4.1 


< 10.5 


KKt] 


6.90 ±0.05* 








8 / 6.6 


< 7.4 


< 3.5 




9.42 ±0.12 


4.56 ± 0.05 


1 09+0.11 


1.82 


07 C-+21.3 

■>J_20.8 


< 65.4 


< 27.4 


KKAtt 


10.38 ±0.12 


4.64 ± 0.06 


'^■•^°-0.07 




n Q + 12.6 
^^•■^-12. 2 


< 20.6 


< 7.7 


KsK3n 


11.65 ±0.13 


4.82 ± 0.05 


^- 'O-0.08 


1.78 


12.9l7;g 


< 23.9 


< 11.4 



*Submode-decay branching fractions are included in efficiencies. 

The modes KKtt^ and KsKtt are actually specific final states of KKtt. There- 
fore, we also determined the product branching fraction for the decay ip{2S) — > 
jrjc{2S),rjc{2S) — > KKn. Figure 1^.371 shows the fit result of the combined final 
states. The results are also listed in Table 13 .Ti 



The statistical significance, defined as \/—2A\nC, where is C is the likelihood and 
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the difference is computed between the standard fit with both signal and background 
and an ahernative background-only fit. Excluding the r]c{2S) —>■ Trvrr/' mode, the 
product of branching fractions B{iIj{2S) — > 7 ric{2S j) x B{ric{2S) — > X) {Bi x B2 in 
Table I3.7P was obtained from 

B{^[j{2S) ^ 7 Vc{2S)) X B{r],{2S) ^ X) = (3.4) 

where Ngig is the number of signal events (yields) determined by fits for non-rj modes 
and Feldman-Cousins table for modes contain r^'s, e is the efficiency, N^(^2S) is the total 
number of ^'(25') decays. For rjc{2S) mrr]' mode additional factors of B{ri' rnir]) 
and B{ri 77) should be considered and the corresponding equation becomes 

BmS) 7 Vc{2S)) X B{7^,{2S) ^ TTvrV) = — ^ ^^^^ ^. (3.5) 

When calculating the upper limits of the branching fraction product we used N^(2S) = 
25.9 X 10^ [76j. The results of the upper limits of branching fraction products at 90% 
confidence level are also listed in Table 13.71 

No evidence for the decay ip{2S) — > 7?7c(25') is observed in any of the these modes. 
The only mode in which an excess above background is present with a statistical 
significance greater than 2a is the An mode. We have investigated other aspects of 
the events in the signal region to determine if the excess constitutes evidence of a 
signal or is just a statistical fluctuation. 

Fi gure [3.381 shows the cos 6 distribution of the photon for the An mode. The points 
are determined by the measured photon energy yields after requiring the photon 
candidate to be within a given cos 9 region. This distribution is fit to the expected 
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1 + 008^6* dependence, giving a x^/dof of "i.l jl^ which corresponds to a confidence 
level of only 16%. 

Another test was to fit the distribution in the constrained photon energy, the value 
returned by the 4C kinematic fit to the ^^{28) ■yric{2S), ric{2S) An hypothesis. 
If the events in the excess are true signal, we would expect the significance to be 
increased. In fact it is reduced to only 2.6cr. Figure 13.391 displays the constrained 
photon energy for the 47r mode. 

These tests suggest that the excess of events in the 7]c{2S) signal region for the Att 
mode is caused by an upward fiuctuation of the background rather than true signal. 

Our conclusion is that there is no statistically significant excess attributable to 
i/^{2S) — > ri(.{2S) in any of the 11 modes considered. It is reasonable to ask whether 
the sum over all exclusive modes shows evidence of an excess, even though the in- 
terpretation of such an excess would be complicated. Figure 13.401 (top) shows the 
measured photon energy distribution summed over all 11 modes (points). It can 
be compared with a "background-only" hypothesis constructed as the sum of the 
distributions from mode-by-mode fits to the MC-determined background histograms 
with no signal component included (solid histogram). A small excess is visible in the 
signal region, but it is clearly not statistically significant. 

Figure 13.401 (top) shows a clear discrepancy between the fitted background and 
the observed yield in the lowest photon-energy bins. As was discussed previously, 
this region has a significant contribution from splitoff showers. In our J/ip study 
(Figure [3.211 through 13.241) . we tested the reliability of the background simulation in 
events that are similar to ric{2S) signal events. In addition to parameterizing the 
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Figure 3.38: The cos^ distribution of the photon candidate for the 47r mode. The 
points are the fit yields after requiring the photon candidate to lie with a particular 
cos 6 bin. The line is the result of fitting the points to 1 + cos^ 6 with only the 
normalization left free. 
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MINUIT Likelihood Fit to Plot 
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Figure 3.39: Constrained photon energy for the final state 'jAtt in the r]c{2S) signal 
region. The points are from the 25.9 M ^^{23) data sample. The dotted line is the 
Breit-Wigner convoluted with the Crystal Ball resolution signal shape. The dashed 
line is the background histogram determined from the 10 times luminosity generic 
ip{2S) and 5 times luminosity continuum MC samples. The solid line is the sum of 
the signal and background. 
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background with a single MC-predicted distribution, we considered the alternative 
of fitting to two separated background components: splitoff showers and showers 
from all other sources. While there was some improvement in fit quality for the 
fits with two independent background parameters compared to fits with a single 
background parameter, it was not statistically compelling and we adopted fits with 
one background parameter for our standard procedure. The discrepancy observed 
in Figure 13.401 (top) may demonstrate that for the full set of 11 signal modes the 
standard MC does not provide a satisfactory description of the background. This led 
us to reconsider whether fitting with two background parameters would give better 
agreement. 

Figure 13.401 (bottom) shows the same data points with the summed histogram 
obtained from background-only fits with two parameters, one controlling the con- 
tribution of splitoff showers and one for other showers. The separate background 
components were constrained to be positive. Because of insufficient statistics to 
determine a second background parameter, the treatment of the rj modes was left 
unmodified. The sum of the background distributions for the modified fits in Fig- 
ure 13.401 (bottom) shows clear improvement in the agreement with the data, sug- 
gesting a deficiency in the MC that is at least partially addressed by the fits with 
two background parameters. Our conclusion that there is no statistically significant 
excess attributable to ric{2S) stands, but this exercise provides a reminder of possible 
systematic sensitivity in our results to the treatment of the background. 

From Figure 13.401 because of the low statistics, we cannot find an explanation 
other than statistical fluctuation for the possible discrepancy between data and MC 
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Figure 3.40: Measured photon energy for all 11 modes in the rjc{2S) signal region. 
The points are from the 25.9 M ip{2S) data sample. The histogram is the 10 times 
luminosity generic ^^{23) and 5 times luminosity continuum MC samples, with the 
contribution of each of the 11 final states determined by its individual 1-parameter 
(top) and 2-paramter (bottom) background-only fits. The red dashed histogram 
shows the events with transition photon candidates tagged as splitoff showers in 
MC. Details about the fits can be found in text. 
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around 70 MeV. 

3.7 Systematic Uncertainties 

While no evidence for the ipi^S) 'yric{2S) transition was observed in any of the 11 
modes we studied, systematic uncertainties were evaluated and incorporated into the 
branching fraction upper limits. 

The systematic uncertainties due to various cuts were studied in the following 
environments. Some systematic uncertainties have been performed in other CLEO 
analyses, and the others can be estimated by studying decays of ip{2S) lXc2 
or ip{2S) — >• jr]c{2S) fits. We used ip{2S) 7Xc2, Xc2 — ^ X decays to study the 
selection criteria, although sufficient statistics were available for only six modes of Xc2 
decays: An, Gtt, KKmr, KK-KTrn^, KKAtt and KsKStt. For some other modes, we 
estimated systematic uncertainties by relating to modes with similar combinations of 
final state particles, e.g., for the KKtt^ mode we applied the corresponding systematic 
uncertainties from the Xc2 KKnmr^ studies. We studied the t]c{2S) signal region 
to assess uncertainties in the ric{2S) signal fitting method. 

3.7.1 Systematic Uncertainties That Apply to All Studied Modes 

Some of the systematic uncertainties were applied to all modes. The uncertainty on 
the total number of 'ip{2S) was 2% [76j. The systematic errors on the trigger, which 
arise from uncertainties in the efficiency of trigger decisions for different final states, 
were taken to be 1% for all modes pTj The systematic error in the efficiency 
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for detecting the transition photon was 2% [911 122] • The results of mode dependent 
systematic uncertainties that were determined by studying ip(2S) 7Xc2 decays are 
hsted in Table [331 while the results determined by studying the ric{2S) signal region 
are listed in Table I3.9[ 

Global Event Selection Systematics 

We conservatively take the respective systematic uncertainties associated with x^/dof 
from the full event vertex fit and the x^/dof from the full event 4C fit to be the 
difference between the results obtained with and without those cuts applied. An 
uncertainty of 1.3% was assigned to the x^/dof from the full event vertex fit for all 
modes. For modes without (with) a Tc^,ri — > 77 decay we assigned a 2.2% (4.0%) 
uncertainty due to the x^/dof from the full event 4C fit. The per mode results can 
be found in Table 13.81 

Resolution Function Systematics 

Resolution function systematics were studied by individually adjusting the Crystal 
Ball parameters a and a by one standard deviation and finding the changes in the 
results of branching fractions from ip{2S) — > 7Xc2 decays. We also performed the 
same procedures for the '4'{2S) 7?7c(2S') decay candidates and obtained consistent 
but smaller results. We assigned 0.8% and 0.6% systematic uncertainties due to the 
uncertainties in a and a, respectively. The per mode results from the Xc2 study are 
listed in Table 13.81 



142 



Signal Region Systematics 

We studied the effect of the lower boundary of the signal region through studying 
ip{2S) —>■ 7Xc2 decays by increasing the lower boundary of the signal region from 
90 MeV to 104-110 MeV, which is about 18 MeV below the mean of the Xc2- This 
amount is the same as the difference between the lower boundary in the ric{2S) signal 
region and the mean of the ric{2S). We assigned a systematic uncertainty of 3.2% to 
all modes. The per mode results can be found in Table 13.81 

The systematic uncertainties due to the high side boundary of the signal region 
were studied directly in the ric{2S) signal region. The changes of the upper limits of 
branching fraction products were determined by varying the upper boundary of the 
E^ signal region by ±8 MeV. The results can be found in Table 13.91 

Hadronic Invariant Mass AM Systematics 

The systematics due to the hadronic invariant mass AM cut were estimated by re- 
moving the cut from the ip{2S) 7?7c(2S') decay search and observing the effect it 
has on the upper limits of branching fraction products. This actually studied the 
change in the background. The results are listed in Table 13.91 

We also studied the ip{2S) IXci decays by adding a requirement of AM < 
175 MeV to simulate the effect of adding a cut of AM at the value 30 MeV below to 
the lower boundary of the r]c{2S) signal region. Results are shown in Table [3T8l and 
the sensitivity is found to be smaller than in the ric{2S) signal region. 
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Background Shape Systematics 

The systematic uncertainties associated with the background shape were estimated 
using the r]c{2S) signal region. In non-r] decay modes, the uncertainties were obtained 
by repeating the signal fits but replacing the histogram background with linear or 
constant functions. For t] modes, in which we do not use fits for upper limit estimates, 
we use the same number of observed events but take the background in the signal 
region down by la, where a is the square root of the number of data events in the 
sideband region. The difference between this and the nominal result was assigned as 
the systematic uncertainty. The results can be found in Table 13.91 

Table 3.8: Mode dependent systematic uncertainties determined from iIj{2S) — > 7Xc2 
decays. The uncertainties are due to the foUowing cuts: (A) x^/dof from full event 
vertex fit; (B) /doi from full event 4C fit; (C) resolution function; (D) lower 
boundary of signal region ranging from = (104, 110) MeV; (E) hadronic 
invariant mass AM. The relative uncertainties are listed in percent. 



Mode 


A 


B 


C 


D 


E 


a 


a 


~ 110 MeV 


AM < 175 MeV 


An 


0.8 


0.5 


0.5 


0.4 


3.2 


0.05 


Gtt 


1.3 


2.2 


0.6 


0.5 


3.0 


0.03 


KKtttt 


1.1 


0.6 


0.5 


0.4 


3.0 


0.2 




~ 


4.0 


0.2 


0.3 


0.7 


0.3 


KKAti 


0.5 


0.2 


0.7 


0.6 


1.8 


0.1 




1.1 


1.1 


0.8 


0.5 


2.5 


0.04 
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Table 3.9: Mode dependent systematic uncertainties determined from studying the 
r)c{2S) signal region. The labels A, B, and C refer to studies of the (A) high side 
boundary of E^^ signal region; (B) hadronic invariant mass AM; (C) background 
shape. The relative uncertainties are listed in percent. 



MnHp 


A 


B 


C 


Att 

T: /I 


11.2 


2 5 




67r 


11.7 


16.6 


20.4 


KKtttt 


13.4 


7.7 


32.6 


KKn^ 


11.0 


26.8 


36.5 


KsKn 


9.1 


1.9 


22.2 


KKn 


13.0 


15.2 


16.9 


TTTrrj 




6.1 


46.5 


inrrj' 




8.9 


24.5 


KKr) 




9.2 


29.8 


KKtttitt^ 


24.3 


22.7 


15.2 


KKAn 


10.5 


0.8 


6.8 


KsK3n 


4.6 


13.2 


17.8 
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3.7.2 Systematic Uncertainties That Apply Only to Specific Modes 

Other systematic uncertainties apply to only specific modes. Table 13.101 shows a 
summary of these uncertainties. 

Tracking and PID Systematics 

For tracking uncertainties we applied 0.3% for pions and 0.6% for kaons based on 
studies of hadronic D decays, which used the full 818 pb~^ ^/'(3770) data and standard 
DTag track quality cuts [93j . 

Based on the study of the first 281 pb"^ ^(3770) data [94], the pion (kaon) PID 
efficiency determined from MC was 0.55% (1.06%) higher than determined from data. 
Therefore, we decreased the MC-determined efficiencies by these amounts per track. 
We conservatively assign 0.2% (0.3%) uncertainty to the pion (kaon) PID. 

The uncertainties due to the same type of tracks are combined with full correla- 
tion. The uncertainties due to different type of tracks are considered uncorrelated, 
since they are determined through independent studies using specially selected event 
samples that have minimal overlap with our sample. 

7r° and 77 — 77 Finding Systematics 

The uncertainty for 7r° finding was determined by using the results of the vr*^ finding 
study which used the full 818 pb^^ ■ip{3770) data sample [95] . 

Our 71^ selection criteria are slightly different from the nominal criteria used in 
that study since we use the CCFIX. The CCFIX makes two types of changes to 
standard photon reconstruction. It increases the photon resolution in MC by 20% 
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Table 3.10: Systematic uncertainties applied only to specific modes. The uncer- 
tainties are due to the following selection criteria: (A) pion and kaon track finding; 
(B) pion and kaon PID; (C) 7r° and rj finding; (D) mass ranges for 77 — *■ Stttttt^ 
and 77' — 7r7r77(77) decays; (E) Ks finding; (F) Rec M(7r7r) for J/ijj rejection; (G) 
M{X -tttt) for J/tp rejection; (H) cosO^^^; (I) d^^trk', (J) Rec M(7r° or 7] child); (K) 
Rec M{r]) for J/ip rejection. The relative uncertainties are listed in percent. 
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0.4 
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2.3 


1.0 


7r7rr/(7r7r7r°) 


1.2 


0.8 


3.5 


UA 








2.4 






0.8 


TTTrr)' 


1.2 


0.8 


3.5 


7.6 


















1.2 


0.6 
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1.8 


1.0 


3.5 
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1.0 
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2.2 
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KK4n 


2.4 


1.4 








1.3 
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KsKSn 
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0.9 
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and increases the reported errors of the photon uncertainty in data and MC by 20% 
[77] . Using the result in Ref . [95] with a loose tt^ pull-mass cut of ±6 is then similar 
to using the CCFIX. Defining Ae = edata/cMC — 1; the result for Ae with a ±6 7r° pull 
mass and E9/E25 OK cuts applied to the 7r° children is (—3.5 ±2.5)%. Interestingly, 
Ae = (—3.5 ± 3.3)% with a ±6 tt^ pull mass and no E9/E25 OK cuts applied to the 
7r° children, suggesting the E9/E25 OK cuts do not effect finding. 

The uncertainty due to ?7 — 77 finding was determined based on the studies of 
il'{2S) rjJ/ip decays in the 27.4 M ip{2S) data sample [96]. This study found 
Ae = (-6.5 ± 1.3)% without the CCFIX but applying the E9/E25 OK cuts, Ae = 
(-5.6±1.3)% without the CCFIX and the E9/E25 OK cuts, Ae = (-3.5±1.3)% with 
the CCFIX and without the E9/E25 OK cuts. This suggests that the MC-determined 
efficiency overestimates the r/ — > 77 finding efficiency by 3.5% and that the E9/E25 
OK cut does not effect 77 — >• 77 finding. 

Based on the interpretations of the studies above, we decreased the MC-determined 
efficiencies by 3.5% and conservatively assigned systematic uncertainties of 3.5% for 
both 7r° and 77 ^ 77 finding. The systematic uncertainty for 77 tttttt^ decays was 
determined by applying the assigned vr and vr^ systematic uncertainties. Similarly, 
the systematic uncertainty for 77' 7i7i'r]{'~f'~f) decays were determined by applying the 
assigned tt and 7/ — >• 77 systematic uncertainties. 

Ks Finding Systematics 

The systematic uncertainty associated with Ks finding has been studied in several 
CLEO analyses [971 [981 EHl EH [73] . We applied an uncertainty of 2% based on the 
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analyses which apply selection criteria similar to ours [88| [99]. 
r] TTTTTT^ and 7]' Selection Systematics 

The uncertainties due to the mass ranges used to select rj —>■ tttttt^ and r/' —>■ -k-kt] 
decays were determined by changing the corresponding rj or rj' mass ranges to ranges in 
which the inefficiencies double with respect to our nominal range. These inefficiencies 
correspond to the mass ranges of \Mt^t^t^o — < 8.5 MeV for rj -k-r-k^ decays and 
IM^^tt" ~ Mr]'\ < 8-5 MeV for the t]' imr] mode. 

TTTT Recoil Mass Systematics 

The uncertainties associated with J ftp rejection based on the tctt recoil mass criteria 
were determined by studying 4'{2S) 7Xc2, Xc2 X decays with the cuts removed 
and observing the change of the branching fractions. 

M{X — tttt) Systematics 

The uncertainties due to the cut on the hadronic invariant mass determined while 
excluding a pion pair M{X — tttt) for J/ip rejection were determined based on decays 
of ip{2S) 7Xc2,Xc2 — ^ by removing the cuts and observing the change of the 
branching fractions. 

cos^^^TT cind d^^trk Systematics 

The systematics due to the cuts on the angle between the momenta of the transition 
photon and the closest pion (cos^^^tt), and the distance between the transition photon 
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and the nearest track (dj^trk) were determined by removing the cut and observing the 
change of the final resuhs. For the 4tt and KKittt modes we studied ip{2S) —>■ '~fXc2, 
Xc2 —>■ X decays. For KKtt^ and TnTrjlmnT^) modes, we apphed the resuhs from the 
ip{2S) 7Xc2, Xc2 X decay studies for the KKnn and 47r modes, respectively. 

Photon Recoil Mass Systematics 

The uncertainties associated with the photon recoil mass cuts, which rejects ip{2S) — > 
IXcj events where the transition photon was used as a child in a reconstructed vr^ 
or 7] ^ 77 decay (Rec M(7r° or r] child)), were determined based on the decay of 
i^i^S) 7Xc25 Xc2 KKnnn^ by applying the cut to this mode and observing 
the change of the branching fraction. We applied this uncertainty to the KKtt^, 
7r7r?7(77), and KKri{'~f'~f) modes. 

T] Recoil Mass Systematics 

The systematic uncertainties due to the rj recoil mass cuts (Rec M{ri)), which were 
used for J ftp rejection, were determined by changing the recoil mass range to double 
the inefficiency with respect to the nominal cut. These inefficiencies correspond to 
mass range |RecM(?7) — M{J/iIj) \ > 85 MeV for the 717177(77) mode and |RecM(?7) — 
M{J/^)\ > 35 MeV for the nnr]{nnn^) mode. 

3.7.3 Summary 

Consideration of the full set of systematic uncertainties listed in Tables 13.81 through 
I3.10[ as well as those due to the trigger, the number of ipi'^S), and transition photon 
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finding, reveals that the majority of the overall systematic uncertainty is due to the 
treatment of the background. These factors are listed in Column C of Table [3^ The 
E.~f signal region (Column A of Table 13^ also shows significant uncertainty associated 
with the high limit. However, this sensitivity is due to the fact that changing 
the upper boundary of the E^^ signal region changes the background shape because 
of the low statistics in the E^ = (66, 94) MeV region. Therefore, the results listed in 
Column A of Table [3l9] are mostly redundant with the background shape uncertainties 
(Column C of Table \3M . Table [3?TT] lists the total systematic uncertainty applied 
to each mode. In general, the individual uncertainties are treated as uncorrelated 
and are combined in quadrature to obtain the overall systematic uncertainties in our 
product branching fraction measurements. 
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Table 3.11: Systematic uncertainties of all modes. 



Mode 


Syst Err (%) 




14.3 




29.1 


KKlTTT 


36.5 




47.2 


KsKn 


24.7 


KKtt 


26.9 


TTTTT) 


48.0 


TTTTrj' 


28.1 


KKq 


32.1 


KKTTTnr^ 


37.4 


KKAn 


14.0 


KsKSn 


23.4 
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CHAPTER 4 



Conclusions and Discussions 

4.1 Summary of Results 

We do not observe the transition ip{2S) jric{2S) in any of the decay modes studied. 
Table HIT] hsts the results of the product branching fractions B{i(j{2S) —>■ 7 rjci^S)) x 
B{ric{2S) — * X) using the nominal value of the rjci^S) full width (V{rjc{2S)) = 14 ± 
7 MeV [22] )■ Even though we may have an excess of events in the 47r mode, we did 
not find any strong evidence for the iIj{2S) 7 rjc{2S) decay. 

The effect of the ric{2S) full width uncertainty was determined as follows. Separate 
signal MC samples with r{ric{2S)) = 7 and 21 MeV were generated in the same 
manner as the nominal MC sample described in Section [3.2.31 The measured photon 
energy distributions were fit in the same manner as the nominal procedure, but with 
the resolution functions determined from these MC samples and the full width of the 
signal shape adjusted according to the full width being investigated. The results using 
r(?7c(2S')) = 7 and 21 MeV, and the linear extrapolation of the product branching 
fraction as a function of r(?7c(2S')), are listed in Table HTTl 

The branching fractions of ric{2S) X may be estimated from ric{lS) — > X 
branching fractions and the ratio of the full width of the t]c{2S) and t]c{1S) if the 
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Table 4.1: Summary for tp{2S) — > ^ric{2S), r]c{2S) X product branching fraction 
results using r{r]c{2S)) = 14 McV. The product Bi x B2 is defined as B{i;{2S) 
7 r]c{2S)) X B{rjc{2S) — > X) and are at the 10~^ level. The first three columns only 
include statistical uncertainties, while the last column includes both statistical and 
systematic uncertainties. All upper limits are at 90% confidence level. 



Mode 


N ■ 
(90% C.L.) 


Corrected 

6(%) 


B, X B2 
(stat only) 


Bi X B2 
(stat±syst) 


47r 


< 64.8 


20.49 ±0.16 


< 12.2 


< 14.0 




< 36.6 


14.22 ±0.14 


< 9.9 


< 12.9 


KKlTTT 


< 35.2 


19.49 ±0.15 


< 7.0 


< 9.5 




< 16.0 


17.76 ±0.14 


< 3.5 


< 5.2 


KsKtt 


< 11.0 


20.40 ±0.15 


< 3.0 


< 3.8 


KKtt 


< 21.9 


7.63 ± 0.04* 


< 11.1 


< 14.1 


TX-Rl] 


< 1.3 


5.68 ±0.05* 


< 2.9 


< 1.3 


nnr]' 


< 4.1 


8.14 ±0.10 


< 11.1 


< 14.2 


KKt] 


< 7.5 


6.47 ±0.05* 


< 4.4 


< 5.8 




< 65.4 


8.74 ±0.11 


< 29.2 


< 40.2 


KKAtt 


< 20.6 


9.93 ±0.11 


< 8.0 


< 9.1 


KsKSn 


< 23.9 


11.39 ±0.13 


< 11.7 


< 14.4 



*Submode-decay branching fractions are included in efficiencies. 
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Table 4.2: Summary of results for ^(25) ^ 7r?c(2S') with r(r/c(2S')) = 7, 21 MeV. 
All upper limits are at 90% confidence level. The efficiencies include PID, 7r° 
finding, and rj finding corrections. The product Bi x B2 is defined as B{iIj{2S) 
7 r]ci2S)) X ^(r?c(2S') X) and are at the 10"^ level. The Bi x B2 results include 
statistical and systematic uncertainties. The y-intercept and slope variables a and 
b are determined by Bi x ^2 < a + 6 * r{r]c{2S)). 



Mode 


r(?7c(25)) = 7 MeV 


r(?7c(25)) = 21 MeV 


a 


b 

(10~^ MeV~^) 




e (%) 


Bi X B2 




e (%) 


Bi X B2 


47r 


< 53.1 


22.06 


< 10.6 


< 77.5 


19.41 


< 17.7 


7.04 


0.505 


Gtt 


< 26.4 


14.71 


< 9.0 


< 49.8 


13.03 


< 19.1 


3.88 


0.727 


KKtttt 


< 25.6 


20.44 


< 6.6 


< 45.7 


17.72 


< 13.6 


3.10 


0.500 


KKtt^ 


< 12.0 


19.15 


< 3.6 


< 19.5 


16.88 


< 6.6 


2.08 


0.217 


KsKtt 


< 9.7 


21.78 


< 3.1 


< 12.4 


19.53 


< 4.4 


2.43 


0.095 


KKtt 


< 17.2 


8.21* 


< 10.3 


< 26.7 


7.31* 


< 17.9 


6.48 


0.542 


TTTrrj 


< 4.3 


6.79* 


< 3.6 


<4.3 


4.97* 


< 4.9 


2.95 


0.095 


TTirr]' 


< 4.1 


9.46 


< 12.3 


<4.1 


6.98 


< 16.6 


10.1 


0.309 


KKr) 


< 7.5 


7.72* 


< 5.0 


< 7.5 


5.68* 


< 6.7 


4.08 


0.127 




< 49.4 


9.47 


< 28.0 


< 83.9 


8.16 


< 55.2 


14.4 


1.95 


KKAn 


< 17.0 


10.50 


< 7.1 


< 24.6 


9.37 


< 11.6 


4.91 


0.317 


KsKStt 


< 20.2 


12.00 


< 11.6 


< 27.4 


10.23 


< 18.4 


8.19 


0.486 



*Submode-decay branching fractions are included in efficiencies. 
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matrix elements governing the decays of the ric{2S) are similar to those of the ?7c(lS'). 
Since the studies of ip{2S) 7 Xc2^ Xc2 ~^ X built confidence in this analysis proce- 
dure, the most likely explanation of the non-observation is that the branching fraction 
assumptions were overly optimistic. Our original objectives for measuring the reso- 
nance properties of the ric{2S) cannot be achieved with this data sample. The upper 
limits of the branching fractions were obtained instead. 

4.2 Comparison with Published Results 

The recent CLEO measurement for B{J/il) 7 ric{lS)) requires a reestimate of 
Bi^ipi^S) 7 ric{2S)), as computed using Equation 1 1.1 II in Section [1.3.4[ Reaveraging 
B{J/tlj 7 r],{lS)), after including the CLEO resuh of B{J/iIj 7 t]c{IS)) = 
(1.98 ±0.09 ±0.30)% [TOO], leads to i3( J/ ^ 7 7/^(15)) = (1.72 ±0.25)%. Using the 
PDG 2008 values for M(r/c(15)), r(J/V'), M{r],{2S)), and T{tlj{2S)) [52j, along with 
this new value for B{J 7 ri(.{lS)), replaces the prediction given in Equation 11.141 
with B{i){2S) 7 r]c{2S)) = (3.9 ± 1.1) x 10-^ which is 50% larger. 

The BaBar collaboration recently measured the branching fraction for t]c{2S) 
KKn to be B{r]^{2S) KKn) = (1.9 ± 0.4 ± 0.5 ± 1.0)% where the first 
two errors are statistical and systematic, respectively, and the last error is from the 
measurement of B{B+ K+r]^{2S) + c.c.) = (3.4 ± 1.8) x 10"^ [M]. Using 
the central value of the BaBar B{ric{2S) —>■ KKtt) result with our 90% confidence 
level upper limit of B{ip{2S) ■yric{2S)) x B{rjc{2S) KKn) leads to an upper 
limit of B{tp{2S) —>■ 'yr]c{2S)) < 7.4 x lO"'*. This limit is roughly twice the revised 



156 



phenomenological prediction given in the last paragraph. 

The BaBar measurement of B{ric{2S) KKtt) can also be used to obtain an im- 
proved estimate of the two-photon partial width of the ric{2S). CLEO [17j previously 
measured the ratio 



and assumed equal branching fractions to KKtt) for ric{lS) and ric(2S) to obtain a 
first estimate of r^^(r/c(25)) = (1.3 ± 0.6) keV. Using BaBar's B{t]^{2S) KKtt), 
the PDG value of T^^{r]^{lS)) = (7.2 ± 2.1) keV, and B{7],{1S) KKtv) = 7.0 ± 
1.2% [52J, the CLEO measurement of R{r]c{2S) / ric{lS)) can be reinterpreted to give 
T^^{7]c{2S)) = (4.8 ± 3.7) keV. 

A recent analysis from the Belle collaboration found r^.y {ric{2S))B{ric{2S) 
47r) < 6.5 eV and T^^{ri^{2S))B{ri^{2S) KKtttt) < 5.0 eV |55j, both at 90% 
confidence level. Using the T^^{ric{2S)) result above, these two results give upper 
limits of B{r],{2S) Atv) < 0.14% and B{r],{2S) KKtttt) < 0.10%. 

Our results for the 4tt and KKtttt modes can be used to determine comparable 
upper limits. Using B{ip{2S) ■yric{2S)) = (3.9 ± 1.1) x lO^'' leads to upper limits 
of B{ri^{2S) Att) < 3.6% and B{ri^{2S) KKtttt) < 2.4%. Our upper limit results 
are an order of magnitude larger than the corresponding upper limits determined 
from the Belle results. 

The branching fraction measurement for ric{2S) — > KKtt and the upper limits 
for ric{2S) Att, KKtttt from the Belle collaboration are all lower than predictions 



R{Vc{2S)/v,{lS)) 



T.,.,{r],{2S)) X B{r],{2S) ^ KKtt) 
r^^{r],{lS)) X BivcilS) ^ KKtt) 



0.18 ±0.05 ±0.02 



(4.1) 
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made by assuming the partial widths of the these decays are the same as for the 
r]c{lS). Using the PDG 2008 values for the ri^{lS) and ri^{2S) full widths, T{r]^{lS)) = 
26.7±3.0 MeV and r(?7c(25')) = 14±7 MeV, and assuming the partial widths for the 
?7c(lS') and r]c{2S) decays are identical, leads to all r]c{lS) branching fractions being 
increased by 1.9±1.0 for r]c{2S). Using the PDG 2008 values for the r]c{lS) branching 
fractions B{r]c{lS) KKn) = (7.0 ± 1.2)%, B{r]c{lS) An) = (1.20 ± 0.30)%, and 
B{r]^{lS) KK-k-k) = (1.5±0.6)%, leads to i3(r/c(25) ^ KKtt) = (13.4±2.3±6.9)%, 
B{ri^{2S) 47r) = (2.3 ± 0.6 ± 1.2)% and B{r]^{2S) ^ KKmr) = (2.9 ± 1.1 ± 
1.5)%, where the last error is from the full width scale factor. The B{ric{2S) 
KKtt) measurement is 7.1 x smaller than the partial width scaling assumption but, 
after considering the uncertainties, they only differ by 1.6 a. The Att and KK-r-r 
upper limit results are 16. 4x and 21x smaller, respectively, but only differ by 1.6 
and 1.5 a. Higher precision measurements of the ric{2S) full width and B{B^ 
K^r]c{2S) + c.c.) are needed in order to determine if the partial width scaling 
assumption is correct for the ?7c(l'S', 2S) states. It should be noted that this may not 
be true since, for example, we see large differences in the decay rates of J/ip and 

^{2S) to pTT. 

Is the experimental measurements for the ric{2S) full width correct? Comparing 
the J/tp and ip{2S) full widths would lead us to think the rjc{2S) full width should be 
larger than the ric{lS). Experimentally, this is not true. Is there something unusual 
with the full widths of the J /ip or %p{2S)l Table compares the full width and e"^e^ 
partial widths for the first two S-wave spin triplet states in the cc and hh systems, 
along with the full widths and 77 partial widths of the ri(.{lS) and ric{2S). The 
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full width ratio of ric{2S)/ric{lS) agrees with T{2S)/T{1S) and not with ip{2S)/ J/ip, 
while the e+e~ and 77 partial widths for all three agree. This seems to suggest that 
the ip{2S) full width is the odd character, possibly related to DD threshold effects. 

Table 4.3: Ratio of widths between 25* and IS* states in quarkonia. The "x" in the 
table denotes the type of partial widths, which is 77 for 7]c{lS), rjc{2S) and e^e~ 
for J/i), ipi2S), T{IS), and T{2S). All values are derived from the PDG 2008 [52]. 





r],{2S)/r],{lS) 


^P{2S)/J/^ 


T{2S)/T{1S) 


r(2^)/r(i5) 


0.52 ±0.27 


3.40 ±0.12 


0.59 ±0.05 


r,(2S)/r,(is) 


0.67 ±0.55 


0.429 ±0.013 


0.457 ±0.010 



4.3 Conclusions and Future Prospects 

In order to experimentally observe decays of ^'(25') —>■ 'yric{2S), our results show that 
a much larger data sample is needed. CESR-c and CLEO-c have ceased operation, 
but the more powerful collider BEPC-II and its detector BES-III are now beginning 
operation in Beijing, China. BES-III is expected to collect billion ip{2S) events in 
one year of running [102j . forty times the CLEO-c ip{2S) sample. Our results are an 
important input for the development of future analyses of 4'{2S) ■yric{2S). 
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APPENDIX A 



Resolution Function Fits 

The resolution function fits of the measured photon energy for the ip{2S) —>■ 7Xc2 and 
ip{2S) jric{2S) decays are given here. A Crystal Ball function [8T|,|82] was used to 
fit the resolution distribution from the signal MC samples as described in Section l3.2.3[ 
The resolution distribution is the histogram of the difference between the measured 
and generator-level photon energy. The RooFit software package [83j was used to 
perform the fits. The resolution function fits for ip{2S) 'yric{2S),ric{2S) — > X 
are shown in Figures lA.ll through IA.61 The resolution function fits for ip{2S) 
lXc2, Xc2 X are shown in Figures I A. 71 through IA.13I Some of the Xc2 fits give 
relatively poor values of x^i reflecting the fact that the Crystal Ball function is an 
imperfect description of the photon energy resolution. The systematic uncertainty 
related to this effect is negligible. 



160 



471 



^200 



§1000 
> 

LU 



800 
600 
400 
200 





RMS= 0.00511210.000032 




Mean = 0.000054 ± 0.000045 








- % 


Entries = 13112 




a= 1-50±r; 
a= 0.00468 ±:~GeV 
N(ii^(2S))= 13049±:iV4' 
mean = 0.00032 ± GeV 




n = 140.0 




X^/dof = 66.61/96 (x^ Fit) 










-. , . , 1 . . . . 1 . . . . 1 , . , . 1 . . . . 1 . . . . 1 . , . . 1 . . . . 1 . . . , 1 , . . . 





-0.05 -0.04 -0.03 -0.02 -0.01 



0.01 0.02 



0.03 0.04 0.05 

E^E,„c(GeV/c=') 



Figure A.l: Resolution function for the decay mode Tp{2S) — > 7r/c(2S'), r/c(25) — > An. 
The points are from the r]c{2S) signal MC and the solid line is the result of the fit 
to the Crystal Ball function. 
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Figure A. 2: Resolution functions for the decay modes ip{2S) j7]c{2S),7]c{2S) 
Gtt (top) and KKnir (bottom). The points are from the r]c{2S) signal MC and 
solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 3: Resolution functions for the decay modes ip{2S) — > '^rjc{2S),rjc{2S) 
KKtt^ (top) and KsKtt (bottom). The points are from the r]c{2S) signal MC and 
the solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 4: Resolution functions for the decay modes ip{2S) — > ^7]c{2S),7]c{2S) 
rnrr] (top) and KKt] (bottom). The points are from the rjc{2S) signal MC and 
solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 5: Resolution functions for the decay modes ip{2S) — > 'yrjc{2S), rjc{2S) 
TTvrry', -q' — > 7r7r7/(77) (top) and KK-KinT^ (bottom). The points are from the r]c{2S) 
signal MC and the solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 6: Resolution functions for the decay modes ip{2S) — > 'yrjc{2S),7]c{2S) 
KK4tt (top) and KsKSn (bottom). The points are from the r]c{2S) signal MC and 
the solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 7: Resolution function for the decay mode 'ip{2S) 7Xc2) Xc2 47r. The 
points are from the Xc2 signal MC and the solid line is the result of the fit to the 
Crystal Ball function. 
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Figure A. 8: Resolution functions for the decay modes ip{2S) 7Xc2,Xc2 — > Svr 
(top) and KKtttt (bottom). The points are from the Xc2 signal MC and the solid 
line is the result of the fit to the Crystal Ball function. 
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Figure A. 9: Resolution functions for the decay modes ip{2S) — > 7Xc2,Xc2 — > KKtt^ 
(top) and KsKtt (bottom). The points are from the Xc2 signal MC and the solid 
line is the result of the fit to the Crystal Ball function. 
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Figure A. 10: Resolution functions for the decay modes 'tp{2S) 7Xc2;Xc2 
7r7rr/(77) (top) and Tnrr]{Tnnr^) (bottom). The points are from the Xc2 signal MC 
and the solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 11: Resolution functions for the decay modes 'tp{2S) 'yXc2,Xc2 

7/(77) (top) ciiid KKrj{'KTi'K^) (bottom). The points are from the Xc2 signal 
MC and the solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 12: Resolution functions for the decay modes Tp{2S) 7Xc2; Xc2 mrr]' , 
r( — > 7r7rr/(77) (top) and KK-Kmr^ (bottom). The points are from the Xc2 signal 
MC and the solid line is the result of the fit to the Crystal Ball function. 
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Figure A. 13: Resolution functions for the decay modes tp{2S) — > 7Xc2, Xc2 KKAn 
(top) and KsK3tt (bottom). The points are from the Xc2 signal MC and the solid 
line is the result of the fit to the Crystal Ball function. 
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APPENDIX B 



Upper Limit Determination Procedures 

Two methods were employed for determining the upper hmits on the yields listed in 
Table 13.71 For modes which do not contain t] decays, we have no apparent signals 
and substantial backgrounds. Therefore, we determined a 90% confidence level upper 
limit by integrating the distribution defined from the nominal yield result up to 90% 
of the area in the physical range (that is with the lower integration limit of zero yield) . 
We used toy MC studies (explained below) to verify that we get consistent results 
with the nominal fit integration. For modes with rj decays, a very limited number 
of events pass our selection criteria, either signal or background. Therefore, we use 
the Feldman and Cousins method [20], which is a standard procedure in high energy 
physics for setting Poisson upper limits in the presence of nonzero backgrounds that 
are separately estimated. 

The distributions of the fitted number of signal events were best represented by 
bifurcated Gaussian distributions. However, to determine the 90% confidence level 
upper limit for each distribution, we find the value of the number of signal events 
at which the area of the distribution between and this value is 90% of the area. 
This procedure is standard in high energy physics and avoids the artificially stringent 
limits that would result from a nonphysical central value. 
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For the modes without rj decays, the measured photon energy were individually 
fit with the signal function and background histogram for all modes, as shown in 
Figures [3.321 through 13.351 For each mode, we used the histogram from the fit result 
as the mean for generating toy MC samples. In a toy MC sample, for each bin 
we use a random number generator to randomly assign the bin content based on 
the Poisson probability distribution with the mean equal to the histogram bin value. 
After all bins of the toy MC sample histogram were filled, a simulated data sample 
was produced. This simulated data sample was then fit with the same signal function 
and background histogram as the real data sample. The signal yield of this fit was 
recorded. For each mode, 10,000 toy MC samples were generated and fit, producing 
a distribution of the 10,000 yields. Figure [BTT] and [B.2I show the yield distributions 
for the 10,000 toy MC samples for all non-r/ modes. Table [RT] compares the toy MC 
results with results from integrating the physical region of the yield from the nominal 
data fit. 

For the modes with t] decays, the following procedure was followed. The signal 
region was defined in the range = [34, 62] MeV, corresponding to ± one full width 
(r(?7c(2S')) = 14 MeV) about the nominal central value. The number of observed 
events is counted. Next, the sideband region defined as = [66, 94] MeV was 
fit with only the background histogram to determine the background normalization 
over the r]c{2S) signal region. Using this background normalization, the number of 
background events in the signal region is estimated. Using the number of observed 
and background events in the signal region, we applied the Feldman and Cousins 
procedure to determine the 90% upper limit on the number of signal events. The 
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Table B.l: Comparison of upper limit methods. The column labeled "Fit Inte- 
gration" is the result from determining the 90% level of the distribution from the 
nominal fit result with only positive area. The column labeled "Toy MC" is the 
result determining the 90% level from the toy MC samples with fit yield > 0. 



Mode 


Fit Integration 


Toy MC 


47r 


< 64.8 


< 65.1 




< 36.6 


< 36.3 


KKlTTT 


< 35.2 


< 34.6 




< 16.0 


< 15.3 


KsKtt 


< 11.0 


< 10.7 


KKttttt:^ 


< 65.4 


< 66.1 


KKAn 


< 20.6 


< 20.8 


KsK37r 


< 23.9 


< 24.5 
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efficiency for tlie signal region, as compared to tlie = [30, 94] MeV range for non-i] 
modes, is tlien used to determine the upper limit on the product branching ratio. 
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Figure B.l: A^'gignai distributions based on the signal and background fit for 10,000 
toy MC samples for the modes Air (top left), KKtttt (top right), KKn^ (bottom 
left), and KsKtt (bottom right). 
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Figure B.2: A'^gignai distributions based on the signal and background fit for 10,000 
toy MC samples for the modes Gvr (top left), KKmnT^ (top right), KKAtt (bottom 
left), and KsKStt (bottom right). 
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